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ABSTRACT 

Context. The primordial nature of the Spite plateau is at odds with the WMAP satellite measurements, implying a primordial Li 
production at least three times higher than observed. It has also been suggested that A(Li) might exhibit a positive correlation with 
metallicity below [Fe/H] — 2.5. Previous samples studied comprised few stars below [Fe/H]=-3. 

Aims. We present VLT-UVES Li abundances of 28 Halo dwarf stars between [Fe/H]=-2.5 and -3.5, ten of which have [Fe/H]< -3. 
Methods. We determined stellar parameters and abundances using four different Teff scales. The direct infrared flux method was 
applied to infrared photometry. Ha wings were fitted with two synthetic grids computed by means of ID LTE atmosphere models, 
assuming two different self-broadening theories. A grid of Ha profiles was finally computed by means of 3D hydrodynamical at- 
mosphere models. The Lii doublet at 670.8 nm has been used to measure A(Li) by means of 3D hydrodynamical NLTE spectral 
syntheses. An analytical fit of A(Li)3D NLTE as a function of equivalent width, T^g, log g, and [Fe/H] has been derived and is made 
available. 

Results. We confirm previous claims that A(Li) does not exhibit a plateau below [Fe/H]=-3. We detect a strong positive correlation 
with [Fe/H] that is insensitive to the choice of T^ff estimator. From a linear fit, we infer a steep slope of about 0.30 dex in A(Li) per 
dex in [Fe/H], which has a significance of 2-3 cr. The slopes derived using the four T^ff estimators are consistent to within lo". A 
significant slope is also detected in the A(Li)-reft- plane, driven mainly by the coolest stars in the sample (reff<6250), which appear 
to be Li-poor. However, when we remove these stars the slope detected in the A(Li) - [Fe/H] plane is not altered significantly. When 
the full sample is considered, the scatter in A(Li) increases by a factor of 2 towards lower metallicities, while the plateau appears very 
thin above [Fe/H]=-2.8. At this metallicity, the plateau lies at (A(Li)3D.NLTE> = 2.199 ± 0.086. 

Conclusions. The meltdown of the Spite plateau below [Fe/H] 3 is established, but its cause is unclear. If the primordial A(Li) 

were that derived from standard BBN, it appears difficult to envision a single depletion phenomenon producing a thin, metallicity in- 
dependent plateau above [Fe/H]=-2.8, and a highly scattered, metallicity dependent distribution below. That no star below [Fe/H]=-3 
lies above the plateau suggests that they formed at plateau level and experienced subsequent depletion. 

Key words, nuclear reactions, nucleosynthesis, abundances - Galaxy: halo - Galaxy: abundances - cosmology: observations - stars: 
Population II 



* Based on observations made with the ESO Very Large Telescope 
at Paranal Observatory, Chile (Programmes 076.A-0463 and 077. D- 
0299). Table 3 is fully available in electronic form at the CDS 
via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http : //cdsweb . u-strasbg . fr/cgi-bin/qcat? J/A+A/ 



1. Introduction 

Spite & Spite (1982a,b) first noted that metal-poor 
(-2.4<[Fe/H]<-1.4), warm (5700 K < reff< 6250 K), dwarf 
stars exhibit a remarkably constant Li abundance, irrespective 
of metallicity and effective temperature, and interpreted this 
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plateau in Li abundance (hereafter the Spite plateau) as being 
representative of the abundance of Li synthesized during the 
primordial hot and dense phase of the Universe (Big Bang, 
Wagoner et al. 1967; see locco et al. 2009 for a review). 
Determining the lithium abundance in unevolved metal-poor 
stars has since developed into an active research topic, because 
of its potential role as a cosmological diagnostic. In the standard 
Big Bang nucleosynthesis (SEEN) scenario, ^Li is formed 
immediately after the Eig Eang, together with 'H, ^H, ^He, and 
^He. is formed first, and is subsequently required as a seed 
to form any heavier element (the so-called "deuterium bottle- 
neck"). The abundance of all the subsequent BBN products thus 
depend on the equilibrium abundance, which is determined 
by the photodissociation reaction ^H(y,'H)'H. As a result, 
all the abundances of BBN products ultimately depend on the 
primordial baryon/photon ratio rjB = nslny (Steigman, 2001), 
and can in principle be employed to constrain this fundamental 
cosmological parameter. 

Following the cosmic microwave background anisotropy 
measurements of WMAP (e.g. Dunkley et al., 2009), tib can be 
inferred from the value of the baryonic density, D.b, i. e. deter- 
mining the primordial abundance of BEN products is no longer 
the only means by which it is estimated. On the other hand, the 
comparison between the two estimates remains of paramount 
importance as a test of the reliability of the BBN theory, of our 
present understanding of the subsequent chemical evolution of 
the elements involved, and, in the case of Li, of our understand- 
ing of stellar atmospheres. 

Among the available BEN products, and ^Li are the most 
reliable rjB indicators. Being never produced in stars, its 
abundance in a low-metallicity environment can be assumed to 
be quite close to the cosmological value. In addition, its sen- 
sitivity to rjB is monotonic and quite strong ( (^H/^H)oc rj'^'^, 
Steigman, 2009). On the other hand, can be effectively mea- 
sured only in high-redshift, low-metallicity damped Lyman a 
(DLAs) or Lyman limit systems, for which the observations are 
so challenging that only seven such high quality measurements 
exist to date, which were all obtained after lOm-class telescopes 
became available (Pettini et al., 2008). The r]B value inferred 
from them is in good agreement with that derived from WMAP 
(Steigman, 2009). 

In contrast, ^Li can be measured with relative ease in the 
photospheres of warm, unevolved stars. The observations are 
typically restricted to dwarfs, at least when one is interested 
in determining the primordial Li abundance, because the frag- 
ile Li nucleus is destroyed by the ^Li(p,a)'*He reaction as soon 
as the temperature reaches 2.6 million K. This implies that giants 
should not be considered, since their deep convective zones mix 
the surface material with layers that exceed this temperature, and 
almost all Li is rapidly destroyed. The ease with which ^Li is de- 
stroyed has always constituted a challenge to existing models of 
convection and diffusion in stellar atmospheres, which predict a 
depletion of at least a factor of four relative to the primordial 
abundance (Michaud et al., 1984). While one could infer that 
some depletion might have occurred, it appeared impossible to 
obtain a constant depletion over such a wide range of effective 
temperatures. The simplest solution was to assume that no deple- 
tion was indeed taking place. This is in marked contrast to the 
solar case, where the photospheric Li abundance is about two 
dex lower than the meteoritic value. 

The original interpretation of the Spite plateau has been chal- 
lenged in many ways in the years since its discovery. Surely the 
most compelling challenge was the independent measurement 
of riB by the WMAP satellite, placing the expected primordial 



Li abundance at A(Li)V=2.65+°°^ (Steigman, 2007), or even 
higher, A(Li)p=2.72±0.05 when updated rates are taken into ac- 
count for the ^He(Q', y)''Li reaction (Cyburt et al., 2008). The 
highest estimate of the Spite plateau does not exceed A(Li)=2.4, 
a more typical value being A(Li)~ 2.2. The discrepancy can in 
principle be eliminated in two ways, by either rejecting the stan- 
dard BEN scenario (for a review see locco et al., 2009), or by 
assuming that some degree of Li depletion has occurred. Two 
main mechanisms could again be invoked. Li could be subject to 
depletion before the currently observed stars are formed (Piau 
et al., 2006), by means of the reprocessing of the primordial 
gas in a first generation of massive, hot stars. This phenomenon 
does not appear to be able to explain the entire WMAP / Spite 
plateau gap, but, removing up to 0.3 dex of the discrepancy could 
considerably reduce the problem. The maximum possible deple- 
tion is nevertheless dependent on the initial mass function and 
lifetime of Pop. Ill stars, as well as on the effectiveness of the 
mixing of their ejecta in the interstellar medium, which are all 
poorly known. Alternatively, Li can be depleted within the stars 
we currently observe, as a consequence of phenomena within the 
envelope, such as difl'usion, gravity waves, rotational mixing, or 
any combination of these. As stated above, the negligible scatter, 
and apparent lack of slope in the Spite plateau are observational 
constraints that models of Li depletion have failed to reproduce. 
This could apparently be achieved by combining diffusion with 
some form of turbulence at the bottom of the atmospheric con- 
vective zone (Richard et al., 2005; Korn et al., 2006, 2007; Piau, 
2008; Lind et al., 2009B). Unfortunately, the effect of turbulence 
is introduced basically as a free parameter, and its tuning is made 
quite difiicult by the subtlety of the efi'ects expected on elements 
other than Li (see sect. 6.4 in Bonifacio et al., 2007). Claims have 
been made (e.g. Asplund et al., 2006) that the lighter *'Li isotope 
has been detected in the atmospheres of dwarf stars displaying 
Spite plateau ^Li abundances. These measurements are very dif- 
ficult and sensitive to subtle details of the analysis (Cayrel et al., 
2007). If the detection of ^Li in EMP dwarf stars were to be con- 
firmed, it would severely undermine any claim of a substantial 
atmospheric depletion of ^Li during the star's lifetime, since the 
^Li is even more easily destroyed than ^Li. 

One additional problem is constituted by repeated claims that 
the Spite plateau might display a tilt towards lower Li abun- 
dances at lower metallicities, on the order of 0.1-0.2 dex in 
A(Li) per dex in [Fe/H] (Ryan et al., 1996, 1999; Boesgaard et 
al., 2005; Asplund et al., 2006), although other studies failed to 
confirm this (e.g. Bonifacio & Molaro, 1997). Roughly below 
[Fe/H] =-2.5, more and more stars appear to exhibit Li abun- 
dances below the plateau level, while the scatter increases. 

The extreme case is possibly represented by the lithium 
abundance upper limit of the hyper-iron-poor subgiant HE 
1327-2326 (Frebel et al., 2008, and references therein), which 
should have A(Li)<0.7 (from ID analysis). The interpretation 
of this result is not straightforward. Even rejecting the inter- 
pretation (Venn & Lambert, 2008) that this star might be a 
chemically-peculiar evolved object, the unusual photospherical 
composition of this star has not yet found a satisfactory explana- 
tion. Were the composition of HE 1327-2326 to be indeed pri- 
mordial, its lack of Li would support the Piau et al. (2006) sug- 
gestion of a pollution by material cycled through massive Pop. 
Ill stars. 

Adopting the Piau et al. (2006) hypothesis, one could then 
envision a scenario in which partial pollution by this astrated 
material induces varying degrees of Li "depletion" in EMP stars 



' A(Li)=Iog[A'(Li)/Af(H)] + 12 
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Table 1. Observations log for the 1 1 new targets. 



Star 


Obs. date 
(UT) 


MJD" 
(UT) 


Exp. time 


vL 


S/N' 


BS 17572-100 


21 Feb 2006 


53787.05628641 


3190 


189 


191 


CS 22188-033 


7 Jul 2006 


53923.31274319 


3035 


14 


107 




6 Till ?no6 


53922.36486180 


4100 

4100 


182 

182 


78 




93 Anr 7006 




3600 


69 


92 


CS 29491-084 


17 May 2006 


53872.31418739 


3600 


-8 


104 


CS 29514-007 


7 Jul 2006 


53923.35103355 


3600 


41 


91 


CS 29516-028 


19 May 2006 


53874.32992778 

53874.37328577 


3600 

3600 


-179 
-179 


62 


CS 30302-145 


23 Apr 2006 


53848.28216715 
53848.30823918 


2160 
2160 


195 
195 


70 


CS 30344-070 


17 May 2006 
27 May 2006 


53872.35925533 
53882.31201609 


2160 
2160 


-141 
-140 


82 


HE 0148-2611 


12 Jul 2006 


53928.37446515 


3600 


-227 


72 


HE 1413-1954 


21 Feb 2006 


53787.27197427 


3600 


-101 


52 



a Modified Julian date of observation start: MJD=JD-2400000.5 
h Rounded to the nearest km/s, barycentric conection applied 
c Near Lii 670.8 nm doublet, if more than one spectrum has been used, S/N 
is measured in the coadded spectrum. 



according to how much this reprocessed gas is available locally 
at the location and time of each stars' formation. A hnear fit to 
EMP stellar Li abundances would then naturally lead to an ex- 
pected trend in A(Li) with [Fe/H], whose slope would appear 
steeper the more the sample is limited to low metallicities. An 
alternative explanation would be to postulate that a Li "over- 
depletion" mechanism operates in the photospheres of the most 
metal-poor stars, a mechanism that would not act uniformly in 
every star of a given metallicity (possibly depending on Tejf or 
rotation speed or both), but would be more efficient at lower 
[Fe/H]. These stars would then begin with a Li abundance cor- 
responding to the Spite plateau, but most of them would then 
develop some degree of Li depletion. This explanation would, 
at the same time, explain the apparent slope at low metallicities 
and the increase in the scatter. It would also explain why, even at 
very low metalUcity, one still finds some stars lying on the Spite 
plateau. A striking example of this is the EMP double-lined bi- 
nary system CS 22876-032 (Gonzalez Hernandez et al., 2008), 
in which, at [Fe/H]=-3.6, the primary lies on the Spite plateau, 
while the secondary has a Li abundance lower by about 0.4 dex. 

2. Observations and data reduction 

Our sample includes 1 1 main-sequence tumofl" and dwarf stars 
selected from various sources (see Table 2), along with the 
sample already presented in Bonifacio et al. (2007). The star 
HE 0148-2611 was previously analyzed (Cohen et al., 2002; 
Carretta et al., 2002) but no Li measurement was ever performed. 
One star (HE 1413-1954) was derived from the Barklem et 
al. (2005) sample. It, again, had no previous Li measurement. 
The remaining targets were drawn from the HK (Beers et al. , 
1985, 1992; Beers , 1999) and Hamburg/ESO (Christlieb et al., 
2008) surveys, and were never studied previously based on high- 
resolution spectra. They were observed by VLT-UVES (Dekker 
et al., 2000) during programmes 076.A-0463(A) (P.I. Lopez, 
HE 1413-1954 and BS 17572-0100) and 077.D-0299(A) (RI. 



Bonifacio, remaining targets). The observation log for the 11 
new targets is in Table 1, where the final signal-to-noise ratio 
(S/N) around the Li 1 670.8 nm doublet is also indicated. For the 
two stars observed during 076.A-0463(A), the observations were 
performed by using VLT-UVES with the DICl dichroic and the 
346nm -I- 580nm setting with a I'.'O slit. These observations thus 
do not contain the 380-480 nm range, but for HE 1413-1954 the 
HERES (Barklem et al., 2005) data were available, which cov- 
ered that wavelength range. For the stars observed during 077. D- 
0299(A), we used DICl with the 390nm -h 580nm setting and 
r.'O slit, thus providing coverage from 360nm to 750imi. All the 
spectra have spectral resolution of R~40 000. The data were re- 
duced using the standard UVES pipehne. In Fig. 1, we show the 
Li I 670.8 nm line region for the 1 1 newly observed stars. 

The data were reduced and analyzed with the same proce- 
dures used in Bonifacio et al. (2007), to which the interested 
reader is referred for details of the analysis and the associated 
uncertainties. An extract of the table listing the employed Fe i 
and Fe n lines, in addition to associated atomic data, equivalent 
widths, and abundances in the 3D scale is available in Table 3. 
The full table is available in the online version. 



3. Atmospliere models and spectrosynthesis 
programs 

3A. 1D LTE models and spectrosynthesis 

Various one-dimensional (ID) local thermodynamical equilib- 
rium (LTE) atmosphere models were employed in the present 
study. F. Castelli's grid of fluxes computed using ATLAS 9 
(Castelli & Kurucz 2003^) was used in the infrared flux temper- 
ature determination (see Sect. 4.2). A second ATLAS 9 model 
grid (Kurucz, 2005; Sbordone et al., 2004; Sbordone, 2005) was 
computed with an ad hoc mixing length parameter in produc- 
ing the Ha-wing profiles used to determine Tetf (see Sect. 4.1). 
It has been shown (Fuhrmann et al., 1993; van't Veer-Menneret 
& Megessier, 1996) that employing in ATLAS a mixing length 
parameter of l/Hp = 0.5 provides the best fit to Balmer fines 
profiles in the Sun, while the value I /Hp = 1.25 generally bet- 
ter reproduces more closely the solar flux, and thus, is usu- 
ally employed in "general purpose" models. To compute these 
profiles, we employed a modified version of R. L. Kurucz's 
code BALMER^, which was capable of handling dififerent line- 
broadening theories. Finally, we employed OSMARCS atmosphere 
models (Gustafsson et al., 1975; Plez et al., 1992; Edvardsson et 
al., 1993; Asplund et al., 1997; Gustafsson et al., 2003) and the 
turbospectrum spectral synthesis code (Alvarez & Plez, 1998) 
to determine Fe i and Fe ii abundances, gravity, microturbulence, 
and Li abundances. Hydrostatic monodimensional LHD models 
(see Caff^au & Ludwig, 2007; CafFau et al., 2007) were used in 
determining the ID NLTE corrections (see Sect. 6.1). 

3.2. 3D hydrodynamical models and spectrosynthesis 

Time-dependent, hydrodynamical 3D stellar atmosphere models 
computed with CO^BOLD (Freytag et al., 2002; Wedemeyer et al., 
2004) as part of the CIFIST model grid (Ludwig et al. 2009B4) 
were employed to produce grids of Ha-wing profiles for T^g es- 



^ grid available at http://wwwuser.oat.ts.astro.it/ 
castelli/ 

' The original version is available online at http: //kurucz. 
harvard . edu/ 

see http://cifist.obspm.fr/ 
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Table 2. Coordinates and optical and infrared photometry for the program stars. 



Star 



H 



K E(B-V) 



BS 16023-046 
BS 17570-063 

BS 17572-100 
CS 22177-009 
CS 22188-033 
CS 22882-027 
CS 22888-031 
CS 22948-093 
CS 22950-173 
CS 22953-037 
CS 22965-054 
CS 22966-011 
CS 29491-084 
CS 29499-060 
CS 29506-007 
CS 29506-090 
CS 29514-007 
CS 29516-028 
CS 29518-020 
CS 29518-043 
CS 29527-015 
CS 30301-024 
CS 30302-145 
CS 30339-069 
CS 30344-070 
CS 31061-032 
HE 0148-2611 
HE 1413-1954 
LP 815^3 



14h00m 
00h20m 

09h 28m 
04h 07m 
00h51m 
OOh 38m 
23h 11m 
21h 50m 
20h 35m 
Olh 25m 
22h 06m 
23h 35m 
22h 28m 
23h 53m 
21h 20m 
21h 30m 
Olh 06m 
22h 25m 
Olh 12m 
Olh 18m 
OOh 29m 
15h 08m 
19h40m 
OOh 30m 
22h 47m 
02h 38m 
Olh 50m 
14h 16m 
20h 38m 



54S.6 
36s. 1 

55S.3 
40S.5 
25S.9 
09S.7 
32S.4 
31S.5 
31S.2 
06S.8 
30S.0 
06S.6 
49S.5 
40S.2 
28S.6 
28S.9 
40S.6 
40S.3 
12S.9 
38S.2 
lOs.5 
29S.7 
52S.2 
15S.9 
23S.2 
43s. 1 
59S.5 
04S.7 
13S.3 



+22° 46' 48" 
+23° 47' 38" 

-05° 21' 36" 
-25° 02' 40" 
-38° 12' 18" 
-31° 47' 54" 
-35° 26' 43" 
-41° 07' 49" 
-15° 53' 30" 
-59° 15' 58" 
-02° 32' 39" 
-30° 22' 53" 
-28° 57' 03" 
-26° 58' 44" 
-20° 46' 24" 
-22° 10' 41" 
-24° 58' 41" 
+05° 37' 40" 
-31° 00' 06" 
-30° 41' 02" 
-19° 10' 07" 
-00° 36' 02" 
-48° 39' 19" 
-35° 56' 51" 
-35° 32' 44" 
+03° 19' 03" 
-25° 57' 02" 
-20° 08' 54" 
-20° 26' 11" 



14.17 
14.51 
12.17 
14.27 
13.20 
15.11 
14.90 
15.18 
14.04 
13.64 
15.10 
14.55 
13.48 
13.03 
14.18 
14.33 
13.97 
15.02 
14.00 
14.57 
14.25 
12.95 
14.46 
14.75 
14.43 
13.90 
14.45 
15.23 
10.91 



13.24 
13.47 
11.28 
13.25 
12.16 



13.02 
13.17 
11.02 
12.96 
11.91 



12.96 
13.07 
10.95 
13.03 
11.90 



0.01801 
0.03949 
0.03727 
0.04407 
0.01315 



13.91 


13.75 


13.65 


0.01417 


14.29 


13.98 


14.00 


0.01576 


12.98 


12.70 


12.66 


0.04551 


12.68 


12.44 


12.46 


0.02796 


13.86 


13.58 


13.45 


0.13321 


13.54 


13.23 


13.27 


0.01391 


12.52 


12.25 


12.20 


0.01367 


12.10 


11.85 


11.86 


0.02027 


13.17 


12.93 


12.87 


0.04547 


13.34 


13.10 


13.07 


0.04547 


12.96 


12.67 


12.66 


0.02375 


13.63 


13.29 


13.15 


0.12816 


13.06 


12.76 


12.74 


0.02241 


13.64 


13.35 


13.37 


0.02030 


13.29 


13.08 


13.05 


0.02213 


11.93 


11.67 


11.67 


0.06527 


13.48 


13.23 


13.26 


0.05343 


13.77 


13.52 


13.45 


0.00904 


13.53 


13.26 


13.27 


0.01305 


12.87 


12.62 


12.61 


0.03727 


13.55 


13.30 


13.30 


0.01362 


14.19 


13.97 


13.89 


0.08681 


9.96 


9.71 


9.65 


0.04514 



timation (see Sect. 4.1, and Ludwig et al., 2009A; Behara et al., 
2009). 

4. Effective temperature 

Effective temperature (T^s) is the most crucial stellar atmo- 
sphere parameter influencing Li abundance determination, Li 
abundances derived from the Lii 670.75nm line being sensitive 
to Teff at the level of about 0.03 dex for each 50 K variation in 
Teff. Unfortunately, a precise determination of stellar effective 
temperatures is generally difficult to achieve. For F/G dwarf and 
subgiant stars such as those studied here, T^ff is routinely esti- 
mated either from photometric calibrations (e.g., Alonso et al., 
2000, 2001) or by fitting the wings of Hor with a grid of synthetic 
profiles of varying T^g. 

Both methods are plagued by specific accuracy issues. 
Photometric calibrations, or the infrared flux method (IRFM), 
are mainly sensitive to the accuracy of the photometry available, 
to the details of the calibration process, and to uncertainties in 
the interstellar reddening estimates. On the other hand, Hor fit- 
ting is mainly sensitive to both the uncertainty in the continuum 
normalization across the broad line wings, and the choice of the 
broadening theory applied in the fine synthesis (see Sect. 4.1). 

In this paper, we considered four temperature estimators: 

- temperatures derived from Hor-wing fitting, using ID at- 
mosphere models and spectrosynthesis, self broadening be- 
ing treated according to Barklem et al. (2000a,b) and Stark 
broadening according to Stehle & Hutcheon (1999) (we will 
henceforth refer to these temperatures as "BA temperatures", 
or the "BA temperature scale"); 

- same as BA, but using the Ali & Griem (1966) self- 
broadening theory (ALI temperatures); 



- temperatures derived from Ha-wing fitting, using 3D model 
atmospheres and spectrosynthesis, Barklem et al. (2000a,b) 
self broadening and Stehle & Hutcheon (1999) Stark broad- 
ening (3D temperatures, for all Hor derived temperatures see 
Sect. 4.1); 

- temperature derived with the infrared flux method (see Sect. 
4.2 as well as Gonzalez Hernandez & Bonifacio, 2009, 
IRFM temperatures). 

4.1. Fitting of ttie Ha wings 

Temperature scales based on Ha-wing fitting are affected by 
both observational and theoretical issues. Most high-resolution 
spectrographs use echelle gratings operating in high orders, 
which exhibit a steep blaze function. The continuum placement 
is thus sensitive to the accuracy with which the shape of the grat- 
ing blaze function can be estimated. Such uncertainties are irrel- 
evant when studying narrow lines observed at high resolution, 
but are important when a broad feature such as Ha is consid- 
ered. More generally, the precision of continuum placement and 
of the determination of the Ha wing shape are affected by noise 
as well as by the possible presence of weak unrecognized fea- 
tures (less of a problem for metal-poor stars). Among these, the 
blaze function shape likely introduces the largest uncertainty. 

On the theoretical side, the uncertainties are due both to 

the atmosphere model structure and to the physics employed in 
the Ha synthesis. Ha-wing self broadening can be treated with 
different theories, most notably those of Ali & Griem (1966), 
Barklem et al. (2000a,b), and Allard et al. (2008). As a gen- 
eral rule, Ali & Griem (1966) theory leads to a significantly 
lower broadening coefficient with respect to the ones derived 
from Barklem et al. (2000a,b) and Allard et al. (2008). A sig- 
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Fig. 1. High-resolution spectra of the Li i 670.79 nm doublet region for 
the 1 1 newly observed stars of the sample. For stars for which multiple 
spectra were available, the coadded spectrum is shown. For the purpose 
of visualization, all the spectra have been shifted to zero radial velocity 
and normalized. T^fi increases from bottom to top, 3D scale T^s and 
A(Li)3DNLTE (see Sects. 4 and 6) are listed for each star. The star CS 
22882-027 shows no detectable Li line, and the 3cr upper limit to A(Li) 
is listed here. 



nificantly higher T^ff is required to reproduce a given observed 
profile when employing the Ali & Griem (1966) theory with re- 
spect to the other theories. With the typical parameters of the 
stars in our sample, and using our fitting procedure, employing 
Ali & Griem (1966) self broadening leads to derived Teff esti- 
mates that are higher by about 150-200 K (a difference of about 
0. 1 dex in Li abundance) with respect to those derived by us- 
ing the Barklem et al. (2000a,b) theory. The theory by AUard 
et al. (2008), on the other hand, leads to T^ff within a few tens 
of K of the Teff estimates obtained when using the Barklem et 
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Fig. 2. Ha red wing profiles in the wavelength range significant for the 
fit. In each panel, red profiles (upper ones) are for reff=5400 K, black 
profiles (middle ones) are for Teff =6000 K, and blue profiles (lower) for 
reff=6600 K. For each temperature, dashed profiles are for log g=3.5, 
solid profiles for log g=4.0, and dotted profiles are for log g=4.5. All 
profiles assume [Fe/H]=-3. Upper panel shows profiles for BA temper- 
atures, lower panel for ALI temperatures. 



al. (2000a,b) theory. We thus restricted ourselves to using the 
self-broadening theory of Barklem et al. (2000a,b) (in BA and 
3D temperatures) and the Ali & Griem (1966) self-broadening 
theory (ALI temperatures). 

The Ha-fitting temperatures exhibit a significant gravity 
sensitivity. Barklem et al. (2002) already reported estimates 
of this sensitivity for relatively metal-poor models (down to 
[Fe/H]=-2). The effect is always in the sense of higher grav- 
ity leading to broader profiles, and appears generally stronger at 
lower metallicities, at lower temperatures, and for the BA pro- 
files compared to the ALI profiles. In Fig. 2, we plot exam- 
ples of profiles for the ALI and BA cases. Profiles are plotted 
for reff=5400, 6000, and 6600 K (higher temperatures generate 
broader profiles). A metallicity of [Fe/H]=-3 is used. For each 
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temperature, we plot the profile for log g=3.5, 4, and 4.5. As can 
be noted, the part of the wing closest to the core appears to be 
more strongly affected than other parts. It is clearly seen that the 
gravity sensitivity of the BA profiles is roughly twice as large as 
in the ALI case. In both the BA and ALI scales, the gravity ef- 
fect becomes quickly negligible as T^ff increases above 6500 K. 
In the most deviant cases ([Fe/H]< -3, reff<6000 K, BA pro- 
files), a difference of 0.5 dex in log g leads to roughly a 200 K 
difference in Teff. 

Thus, the shape of the Ha profile varies in different ways 
when varying gravity and temperature. As a consequence, the 
use of an incorrect value of gravity will always affect the temper- 
ature estimate, but the size of the effect will depend on the details 
of how the actual fitting is performed. To provide some insight 
into what the effect is when employing our specific fitting pro- 
cedure, we fed the fitting program with log g=4, and [Fe/H]=-3 
and -2.5 theoretical Ha profiles, and derived the temperature 
by assuming that log g=3.5, 4.0, and 4.5. When fitting profiles 
of log g=4.0 by means of profiles of log g=3.5, which are nar- 
rower at each temperature, we obtain a higher temperature esti- 
mate than we would if we were to use the proper gravity. The 
opposite effect occurs when using the broader log g=4.5 pro- 
files. In Fig. 3, we plot these temperature differences versus the 
true Teff of the profiles. Differences are computed in the sense 
Areff=reff(log g=3.5 or 4.5) -reff(log g=4.0). Red Hnes with 
filled circles correspond to fits with log g=3.5 profiles, blue ones 
with open diamonds to fits with log g=4.5 profiles. The solid 
lines correspond to [Fe/H]=-3.0 profiles (both fitted and fitting), 
while the dot-dashed line corresponds to [Fe/H]=-2.5. For the 
parameter space covered, and when adopting our fitting proce- 
dure, underestimating the fitting-grid gravity by 0.5 dex leads to 
an overestimate of the temperature by as much as 250 K in the 
BA case, and 200 K in the ALI case. This underestimate reaches 
a maximum around 5200-5300 K, decreasing on both sides, and 
fading away on the hot side, near reff=6500 K. By overestimat- 
ing the fitting-grid gravity, one underestimates T^g by as much 
as 300 K in the BA case and 200 K in the ALI case. The shape 
of the curve is similar, but the point of maximum sensitivity oc- 
curs between 5800 and 6000 K. There is a hint that the effect 
decreases mildly at [Fe/H]= -2.5, although higher metallicities 
have not been explored. 

Since we estimate surface gravity from the Fe i-Fe ii ioniza- 
tion equilibrium, the derived gravity is temperature sensitive, so 
that the two estimations need to be iterated to convergence. As 
a general rule, we stopped iterating when T^ff variations became 
lower than 50 K, which typically required not more than 3 itera- 
tions, starting from an initial guess of log g=4. 

A very mild metallicity sensitivity is also present in the Ha- 
based temperature determination, never surpassing some tens of 
K for a 0.5 dex of variation in [Fe/H]. The actual iteration of the 
temperature determination with the other atmosphere parameters 
was performed differently for the BA and ALI cases on one side, 
and for the 3D case on the other side: 

- In the BA and ALI case, once the gravity and metallicity 
were determined with one temperature estimate, the Ha pro- 
file grid was interpolated to that gravity value, while the near- 
est grid step was chosen in metallicity, without interpolation. 
The small metallicity step of the grid (0.25 dex), as well 
as the very mild sensitivity of Ha to metallicity, made this 
choice sufficiently precise. 

- In the 3D case, the computation of both model atmospheres 
and spectral synthesis is very time consuming, and only the 
atmosphere model grids for log g=4 and log g=4.5 with 
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Fig. 3. Ha theoretical profiles for models with log g=4, [Fe/H]=-3 
(solid lines) and log g=4, [Fe/H]=-2.5 (dot-dashed lines) have been fit- 
ted with the same procedure used for the program stars, but using a log 
g=3.5 grid (red lines with filled circles) and a log g=4.5 grid (blue lines 
with open diamonds). Here we plot the temperature difference (T^ailog 
g=3.5/4.5)-rc(f (log g=4.0) ), against the "real" effective temperature of 
the profile. The upper panel shows BA profiles, the middle panel ALI 
profiles, and lower panel 3D profiles. The gray shaded areas indicates 
the temperature ranges for the program stars in each T^ff scale. 



[Fe/H]=-3.0 were sufficiently extended at the time of the 
analysis. We thus fitted the observed Ha lines to these two 
grids, deriving, for each star, effective temperatures coiTe- 
sponding to the two assumed gravities. We determined grav- 
ity and metallicity, then derived a new Teff estimate by lin- 
early interpolating between the retf(log g=4.0) and re(f(log 
g=4.5) at the estimated gravity. The procedure was then iter- 
ated but, for most stars, the same convergence criterion ap- 
plied to the ID case (AT<50 K) was found to be too strin- 
gent, since the parameters for most stars ended up oscillating 
between two sets corresponding to T^ff estimates that were 
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Table 4. Atmosphere parameters for the program stars using the different temperature estimators. Parameters for CS 22882-027 were derived only 
for the 3D temperature scale. 
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Fig. 4. Effective temperatures for different estimators, plotted against 
BA temperature for the program stars. Top to bottom: ALI, IRFM and 
3D temperatures. The red line represents the one-to-one relation (hence 
the line of BA temperatures). 
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about 60 K apart. This can probably be attributed to the use 
of a more coarse grid in the 3D case. 



Fig. 5. Effective temperatures for different estimators, plotted against 
[Fe/H] (as derived using the temperature in the panel). Top to bottom: 
BA temperatures, then ALI, IRFM and 3D. 
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4.2. IRFM temperature estimation 

Originally introduced by Blackwell & Shallis (1977), and later 
improved by Blackwell et al. (1980), who removed an unneces- 
sary iteration (see Blackwell et al., 1990, and references therein), 
the infrared flux method (IRFM) relies on the ratio of the flux at 
a near-infrared (NIR) wavelength or in a NIR band, to the bolo- 
metric flux. This ratio can also be derived from model atmo- 
spheres, and the efi'ective temperature determined by finding the 
eff'ective temperature of the model that reproduces the observed 
ratio. Gonzalez Hernandez & Bonifacio (2009) presented a new 
implementation of the method making use of 2MASS photom- 
etry (Skrutskie et al., 2006), and also provided a calibration of 
bolometric fluxes with colors (V - J), (V - H), and (V - K^), 
where V is in the Johnson system and the NIR magnitudes are 
in the 2MASS system. We applied the IRFM in exactly the way 
described by Gonzalez Hernandez & Bonifacio (2009), deriving 
the bolometric fluxes as the average of those estimated from the 
three visible-NIR colors. AH magnitudes and colors used in the 
IRFM must be corrected for reddening. To do so, we used the 
reddening maps of Schlegel et al. (1998), corrected as described 
in Bonifacio et al. (2000). All our program stars are sufficiently 
distant that they lie outside the dust layer, so that the full red- 
dening derived from the maps should be applied. The adopted 
reddenings are provided in Table 2. The star CS 22882-027 does 
not appear in the 2MASS catalog, thus we could not derive its 
IRFM temperature. 
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Fig. 6. Difference between A(Li)/ and A(Li), plotted against and 
the Li doublet EW. The two "outUers" are labeled. 



5. Gravity, microturbulence and metallicity 

The FITLINE code was employed to measure the equivalent 
widths of the Fei and Fen lines. Although up to ~120 Fei 
lines were available, only four Fe ii lines were strong enough 
to be used. For each temperature scale, gravity was then derived 
by enforcing Fe i-Fe n ionization equihbrium. For the Hor-based 
scales, gravity was used with metallicity (Fei abundance) to it- 
erate the Teff estimation (see Sect. 4.1). 

For each temeperature scale, microturbulence was deter- 
mined by ensuring that the weak and strong Fe i lines provide the 
same abundance. The final parameters and the derived metaUic- 
ity for each temperature estimator are presented in Table 4, de- 
tailed Fei and Feii abundances are listed in Table 5. Final Teff 
values for the ALI, IRFM, and 3D temperature scale are plotted 
against the BA scale in Fig. 4, and against the respective value 
of [Fe/H]inFig.5. 

6. Lithium abundance determination 

We determined Li equivalent widths in a similar fashion to 
Bonifacio et al. (2007). Synthetic line profiles were fitted to 
the observed profile, and the equivalent width (EW) determined 
from the fitted synthetic profile. The EW errors Usted in Table 6 
were obtained by means of Monte Carlo simulations, in which 
Poisson noise was added to a synthetic spectrum to ensure that 
it had the same S/N as the observed spectrum. The Li abundance 
was determined by iteratively computing synthetic spectra of the 
Li doublet until the synthetic EW matched the observed EW to 
better than 1%. The adopted atomic data were unchanged with 
respect to Bonifacio et al. (2007), and took account of hyper- 
fine structure and isotopic components (a solar Li isotopic ratio 
was assumed). We henceforth refer to these abundances as "ID 
Li abundances" since Li abundances were derived using ID at- 
mosphere model and spectrosynthesis codes. One should avoid 
confusion with the 3D temperature scale, which indicates only 



that 3D effects have been taken into account in determining the 
Hor-wing fitting temperature. 

In addition, we determined, for the 3D temperature scale 
only, what we refer to as "3D NLTE Li abundances". As de- 
scribed in Sect. 6.1, a grid of time-dependent 3D NLTE spec- 
trosyntheses have been produced for the Li i 670.8 nm doublet, 
and used to independently determine Li abundances from the 
measured EW. 

The uncertainties in the Li abundance measurements were 
largely dominated by the uncertainty in the temperature estima- 
tion. For further details, the reader is referred to Bonifacio et al. 
(2007). For the purpose of our analysis, a constant uncertainty 
of o"A(Li) = 0.09 was assumed. 

6.1. 3D and NLTE corrections 

We originally planned to determine the eff'ects of both atmo- 
sphere hydrodynamics and any departure from LTE in a con- 
sistent manner, and thus computed a set of time-dependent 3D 
NLTE syntheses of the Li doublet over a grid of suitable 3D 
models, to construct a set of curves of growth (COG) for the 
doublet EW. Details of the computation of the 3D NLTE lithium 
doublet synthesis are covered in Appendix A. 

The model parameters covered by the COG grid are listed 
in Table 7. Once the COG grid was computed, we decided to 
also derive 3D NLTE lithium abundances directly, by identifying 
the EW-to-abundance relation that most closely fitted the com- 
puted values, and applying it to our observed EW. This was ac- 
complished by either interpolating in the T^g, log g, [Fe/H], and 
EW grid (and possibly extrapolating out of it), or by determin- 
ing a best fitting analytical function in the form A(Li)=f(reff, log 
g, [Fe/H], EW) and applying it to the observed parameters and 
hthium doublet EW. We pursued both of these approaches. 

We found the functional-fit method to be the superior of the 
two, both because of its higher accuracy and for greater ease 
of use. Its functional-form approach condenses the 3D NLTE 
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star 


Fei 


cr 


Fell 


cr 


Fei 


cr 


Fen 


cr 


Fei 


cr 


Fe II 


cr 


Fei 


cr 


Fen 


cr 






BA 






ALI 






IRFM 






3D 




BS 16023-046 


4.53 


0.108 


4.50 


0.135 


4.66 


0.113 


4.64 


0.135 


4.68 


0.115 


4.64 


0.135 


4.56 


0.108 


4.58 


0.136 


BS 17570-063 


4.45 


0.136 


4.51 


0.062 


4.71 


0.156 


4.66 


0.062 


4.64 


0.151 


4.62 


0.061 


4.58 


0.144 


4.61 


0.062 


BS 17572-100 


4.75 


0.124 


4.71 


0.091 


4.88 


0.158 


4.89 


0.085 


4.98 


0.172 


5.02 


0.086 


4.78 


0.152 


4.84 


0.086 


CS 22177-009 


4.33 


0.101 


4.37 


0.123 


4.51 


0.121 


4.54 


0.126 


4.54 


0.114 


4.55 


0.127 


4.42 


0.116 


4.46 


0.123 


CS 22188-033 


4.47 


0.104 


4.44 


0.030 


4.65 


0.123 


4.66 


0.032 


4.52 


0.097 


4.49 


0.032 


4.53 


0.110 


4.57 


0.027 


CS 22888-031 


4.03 


0.148 


4.04 


0.125 


4.32 


0.164 


4.30 


0.132 


4.43 


0.173 


4.35 


0.137 


4.17 


0.154 


4.21 


0.126 


CS 22948-093 


4.19 


0.133 


4.20 


0.144 


4.35 


0.132 


4.31 


0.146 


4.32 


0.129 


4.38 


0.145 


4.26 


0.134 


4.26 


0.146 


CS 22950-173 


4.72 


0.121 


4.76 


0.111 


4.89 


0.126 


4.90 


0.1 10 


4.77 


0.103 


4.77 


0.110 


4.81 


0.124 


4.85 


0.111 


CS 22953-037 


4.59 


0.126 


4.62 


0.124 


4.75 


0.145 


4.74 


0.125 


4.74 


0.136 


4.72 


0.125 


4.66 


0.128 


4.69 


0.125 


CS 22965-054 


4.60 


0.136 


4.62 


0.089 


4.72 


0.141 


4.71 


0.091 


4.71 


0.140 


4.71 


0.091 


4.64 


0.138 


4.67 


0.083 


CS 22966-011 


4.28 


0.106 


4.30 


0.120 


4.54 


0.131 


4.53 


0.121 


4.49 


0.119 


4.48 


0.121 


4.41 


0.118 


4.44 


0.119 


CS 29491-084 


4.46 


0.130 


4.49 


0.119 


4.60 


0.134 


4.58 


0.119 


4.56 


0.133 


4.57 


0.119 


4.53 


0.132 


4.53 


0.119 


CS 29499-060 


4.84 


0.114 


4.82 


0.099 


4.94 


0.127 


4.96 


0.097 


4.94 


0.113 


4.98 


0.100 


4.88 


0.120 


4.89 


0.100 


CS 29506-007 


4.62 


0.116 


4.67 


0.183 


4.80 


0.125 


4.77 


0.186 


4.79 


0.124 


4.77 


0.186 


4.69 


0.118 


4.72 


0.186 


CS 29506-090 


4.67 


0.110 


4.71 


0.109 


4.83 


0.109 


4.82 


0.113 


4.87 


0.122 


4.82 


0.113 


4.73 


0.114 


4.76 


0.109 


CS 29514-007 


4.70 


0.130 


4.66 


0.084 


4.84 


0.136 


4.80 


0.087 


4.71 


0.132 


4.75 


0.078 


4.74 


0.133 


4.74 


0.086 


CS 29516-028 


3.98 


0.135 


3.97 


0.184 


4.31 


0.151 


4.24 


0.183 


4.11 


0.141 


4.10 


0.184 


4.17 


0.150 


4.18 


0.184 


CS 29518-020 


4.64 


0.092 


4.62 


0.138 


4.83 


0.115 


4.83 


0.139 


4.90 


0.125 


4.87 


0.141 


4.71 


0.100 


4.74 


0.139 


CS 29518-043 


4.25 


0.126 


4.28 


0.129 


4.40 


0.128 


4.36 


0.133 


4.34 


0.124 


4.30 


0. 1 33 


4.33 


0.125 


4.31 


0.133 


CS 29527-015 


3.97 


0.143 


3.98 


0.179 


4.13 


0.143 


4.15 


0.179 


4.19 


0.142 


4.19 


0.179 


4.01 


0.145 


4.09 


0.179 


CS 30301-024 


4.79 


0.118 


4.73 


0.106 


4.90 


0.119 


4.92 


0.108 


4.90 


0.119 


4.92 


0.108 


4.81 


0.125 


4.85 


0.102 


CS 30302-145 


4.48 


0.160 


4.48 


0.102 


4.62 


0.160 


4.58 


0.107 


4.62 


0.160 


4.58 


0.107 


4.56 


0.163 


4.53 


0.103 


CS 30339-069 


4.41 


0.188 


4.39 


0.098 


4.57 


0.179 


4.56 


0.101 


4.52 


0.177 


4.55 


0.096 


4.46 


0.189 


4.51 


0.096 


CS 30344-070 


4.48 


0.121 


4.51 


0.104 


4.65 


0.123 


4.61 


0.101 


4.68 


0.125 


4.65 


0.101 


4.58 


0.122 


4.56 


0.102 


CS 31061-032 


4.88 


0.137 


4.88 


0.062 


5.02 


0.132 


4.98 


0.067 


4.92 


0.123 


4.89 


0.062 


4.94 


0.125 


4.93 


0.068 


HE 0148-2611 


4.32 


0.112 


4.31 


0.048 


4.44 


0.118 


4.41 


0.051 


4.43 


0.116 


4.41 


0.051 


4.38 


0.117 


4.36 


0.050 


HE 1413-1954 


4.00 


0.118 


4.01 


0.167 


4.11 


0.118 


4.13 


0.167 


4.28 


0.118 


4.27 


0.168 


4.03 


0.120 


4.09 


0.168 


LP 815^3 


4.62 


0.088 


4.58 


0.051 


4.69 


0.089 


4.71 


0.052 


4.73 


0.091 


4.71 


0.052 


4.70 


0.090 


4.67 


0.052 



Table 6. Lii 670.8 nm EW and errors, and Uthium abundances using the different parameter sets. For the BA, ALI, and IRFM temperature scales, 
we hst ID LTE and ID NLTE A(Li). For the 3D temperature scale, we list ID LTE and NLTE, as well as 3D NLTE A(Li). 



Star 


EW 


error 


A(Li) 


A(Li) 


A(Li) 


A(Li) 


A(Li) 


A(Li) 


A(Li) 


A(Li) 


A(L1) 




pm 


pm 


(a) 


(b) 


(a) 


(b) 


(a) 


(b) 


(a) 


(b) 


(c) 










BA 


ALI 


IRFM 




3D 




BS 16023-046 


1.93 


0.06 


2.145 


2.138 


2.271 


2.257 


2.292 


2.278 


2.193 


2.181 


2.179 


BS 17570-063 


1.76 


0.04 


1.930 


1.928 


2.148 


2.132 


2.091 


2.077 


2.038 


2.025 


2.029 


BS 17572-100 


1.81 


0.04 


2.152 


2.149 


2.232 


2.219 


2.340 


2.329 


2.184 


2.174 


2.166 


CS 22177-009 


2.42 


0.03 


2.153 


2.153 


2.309 


2.293 


2.339 


2.323 


2.224 


2.214 


2.209 


CS 22188-033 


0.78 


0.05 


1.577 


1.577 


1.750 


1.735 


1.665 


1.657 


1.648 


1.636 


1.665 


CS 22888-031 


1.87 


0.04 


1.846 


1.851 


2.104 


2.081 


2.214 


2.199 


1.961 


1.940 


1.976 


CS 22948-093 


1.19 


0.06 


1.935 


1.930 


2.047 


2.034 


2.051 


2.039 


1.988 


1.977 


1.989 


CS 22950-173 


2.11 


0.09 


2.199 


2.193 


2.307 


2.293 


2.212 


2.205 


2.250 


2.238 


2.230 


CS 22953-037 


1.95 


0.03 


2.151 


2.145 


i.ni 


2.258 


2.278 


2.265 


2.211 


2.200 


2.194 


CS 22965-054 


2.21 


0.06 


2.161 


2.164 


2.263 


2.255 


2.277 


2.268 


2.206 


2.203 


2.185 


CS 22966-011 


1.37 


0.05 


1.788 


1.792 


1.983 


1.966 


1.945 


1.929 


1.867 


1.855 


1.869 


CS 29491-084 


1.77 


0.07 


2.080 


2.083 


2.189 


2.180 


2.171 


2.163 


2.142 


2.138 


2.128 


CS 29499-060 


2.07 


0.06 


2.201 


2.196 


2.322 


2.309 


2.311 


2.297 


2.250 


2.239 


2.231 


CS 29506-007 


2.05 


0.04 


2.149 


2.146 


2.275 


2.262 


2.298 


2.285 


2.223 


2.214 


2.205 


CS 29506-090 


1.85 


0.05 


2.102 


2.098 


2.225 


2.211 


2.274 


2.261 


2.156 


2.147 


2.137 


CS 29514-007 


2.33 


0.09 


i.m 


2.209 


2.320 


2.307 


2.231 


2.222 


2.263 


2.254 


2.241 


CS 29516-028 


2.41 


0.12 


1.904 


1.919 


2.159 


2.136 


2.000 


1.990 


2.025 


2.004 


2.035 


CS 29518-020 


2.10 


0.11 


2.052 


2.053 


I.m 


2.195 


2.276 


2.260 


2.109 


2.099 


2.091 


CS 29518-043 


1.72 


0.11 


2.121 


2.115 


2.238 


2.225 


2.204 


2.193 


2.193 


2.182 


2.183 


CS 29527-015 


1.86 


0.06 


2.091 


2.098 


2.188 


2.177 


2.272 


2.259 


2.121 


2.116 


2.113 


CS 30301-024 


1.77 


0.06 


2.143 


2.140 


2.213 


2.200 


2.250 


2.237 


2.157 


2.148 


2.139 


CS 30302-145 


1.54 


0.11 


2.086 


2.077 


2.203 


2.190 


2.215 


2.203 


2.143 


2.131 


2.138 


CS 30339-069 


2.04 


0.11 


2.125 


2.130 


2.223 


2.212 


2.204 


2.194 


2.157 


2.151 


2.140 


CS 30344-070 


1.68 


0.09 


2.064 


2.063 


2.177 


2.166 


2.231 


2.218 


2.135 


2.128 


2.120 


CS 31061-032 


2.10 


0.06 


2.221 


2.213 


2.336 


2.322 


2.250 


2.241 


2.282 


2.269 


2.262 


HE 0148-2611 


1.29 


0.09 


2.000 


1.996 


2.100 


2.089 


2.113 


2.102 


2.065 


2.056 


2.063 


HE 1413-1954 


1 .59 


0.14 


2.035 


2.047 


2.129 


2.123 


2,285 


2.274 


2.079 


2,080 


2.070 


LP 8 15-43 


1.89 


0.06 


2.229 


2.228 


2.302 


1 lyi 


2.334 


2.324 


2.304 


2.296 


2.296 



i ID vameViti'NLTE correction applied. 
c 3D NLTE calculation 



abundance determination into a formula that can be hard-coded 
into any program, eliminating the need to carry over the true 
grid of computed points. Details on the fit calculation, as well 
as the chosen functional form and coefficients, are available in 
Appendix B^. 



' The corresponding IDL functions are also available on-line at 
http : //mygepi . obspm. fr/~ sbordone/fitting .html or by email 
request to the authors. 



As mentioned above, we also performed an interpolation 
over the COG grid. The main problem in producing a suitable 
interpolation lies is that the grid is non-rectangular, which has 
two different causes. First, in the CO^BOLD hydrodynamical mod- 
els Teff is not set a priori, rather, the entropy of the material 
entering through the bottom of the computational box is the 
fixed quantity. The true T^s is determined after snapshot selec- 
tion, and usually varies across an interval of +100 K centered 
on the desired value for the models we employed. As a conse- 
quence, it is impossible to build a grid of CO^BOLD models with 
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Fig. 7. Same as in Fig. 6 but now plotting against log g and [Fe/H]. 



exactly the same temperature but, for example, different metal- 
licity. Secondly, varying the stellar parameters naturally alters 
the relationship between A(Li) and EW, so that the range of EW 
in the COG corresponding to interesting values of A(Li) will 
vary from model to model. 

To simplify the task, we took advantage of the limited sensi- 
tivity of the 670. 8nm Li doublet to both gravity and metallicity. 
Thus, we decided to assume [Fe/H]=-3 throughout the interpo- 
lation, and to avoid interpolating in gravity by always choosing 
the closest value to the derived gravity between log g=4 and log 
g=4.5. This choice was also justified by the limited extension in 
both parameters of our sample. This reduced the problem to in- 
terpolating in an irregularly spaced two-dimensional grid in T^ff 
and EW. Delaunay triangulation^ and quintic polynomial inter- 
polation were then used to derive A(Li). 

Figures 6 and 7 show the difference between A(Li) as de- 
termined by means of the analytical fit {A{Li)f) and by interpo- 
lation (A(Li),), respectively, plotted against relevant quantities, 
on the 3D temperature scale. Most stars show an excellent con- 
cordance between the two methods, but two outliers exist, CS 
29516-028 and CS 22888-031. These two stars have the lowest 
temperatures among all the stars in the sample (these tempera- 
tures are still within the computed grid). However, they also have 
the highest gravity in the sample, which requires extrapolation, 
since the grid has a limiting gravity of log g of 4.5. While the 
functional fit is indeed extrapolated, the simplified interpolation 
assumes log g =4.5 in this instance; the discrepancy between the 
two methods does not however exceed 0.023 dex in A(Li), which 
is negligible for our purpose. All the remaining stars exhibit dis- 
crepancies not exceeding 0.01 dex. 

* Delaunay triangulation is a method of triangulation of a set P of 
points in a plane defined as the triangulation for which no element of 
P lies within the circumcircle of each triangle, except for the triangle 
vertexes. It is often used to model surfaces that are sampled on irreg- 
ular grids (e.g., elevations in geography). The built-in IDL functions 
triangulate and trigrid have been used to produce the triangula- 
tion and the interpolation based on it. 



Fig. 8. The NLTE correction (A(Li)io,A>Lr£ - A(Li) id.lte), computed 
for each star using our model atom (black filled dots) along with the 
Carlsson et al. (1994) values (red open circles), plotted against T^g and 
log g. The 3D temperature scale is assumed. 



A parallel grid of COG was produced using LHD models shar- 
ing the same parameters as the CO''BOLD ones. The ID syntheses 
were produced both including and neglecting NLTE effects, for 
the specific purpose of deriving a grid of NLTE corrections ap- 
plicable to our ID Li abundances. For comparison. Fig. 8 shows 
our ID NLTE corrections (for the 3D temperature scale) versus 
Teff and log g, together with the corresponding values obtained 
by using the Carlsson et al. (1994) NLTE coiTections, while Fig. 
9 shows a similar comparison using the updated calculations by 
Lind et al. (2009). Since Lind et al. (2009) corrections are de- 
fined down to [Fe/H]=-3, in Fig. 9 [Fe/H]=-3 is assumed for all 
stars both in computing our NLTE correction and in computing 
those based on Lind et al. (2009) scale. Trends with effective 
temperature and gravity are extremely similar for our coiTec- 
tions and those of Lind et al. (2009), but a very uniform offset of 
about 0.03 dex is present between the two set of corrections. The 
origin of this offset is probably the different sets of underlying 
atmosphere models. Since the offset is quite uniform across the 
sample, using either set of coiTections is of no consequence on 
the scientific output of the present work. 

7. Results 

We decided to adopt the 3D temperature scale (and its derived 
parameters) together with the 3D NLTE Li abundance set as our 
preferred values, and henceforth, when not otherwise specified, 
we will refer to these. 

7.1. Sensitivity to tine adopted Teff sca/e 

One of the most remarkable results of this work is that, although 
the choice of temperature scale alters the parameters and derived 
Li abundance of the stars, it does not change the general picture 
that emerges. Table 8 provides the results of Kendall's r-test and 
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Fig. 9. The NLTE correction (A(Li)iojv;,j-£ - A(Li) ]d,lte), computed 
for each star using our model atom (black filled dots) along with the 
Lind et al. (2009) values (red open circles), plotted against T^ti and log 
g. The 3D temperature scale is assumed, and [Fe/H]=-3 is imposed for 
all stars. 



the slopes of linear fits to the A(Li)-[Fe/H] and A(Li)-reff re- 
lations. The linear fits were obtained taking into account errors 
in both variables using the fitexy routine (Press et al., 1992). 
The A(Li) error was assumed to be fixed at 0.09 dex, the er- 
ror in [Fe/H] to be given by the Fe i line-to-line scatter for each 
star, and the error in Tgff to have a constant value of 130 K. The 
sample consists of 27 of the 28 stars for which we have Li mea- 
surements, excluding CS 22882-027, for which we only have an 
upper limit to A(Li), as well as CS 21 188-033, where the resid- 
uals from the best-fit regressions are on the order of 3-4cr. 

In statistical terms, for all four temperature scales a non- 
parametric Kendall's r-test indicates that A(Li) correlates with 
[Fe/H] at a very high level of significance (see Table 8). 
Moreover, a linear fit to the A(Li) - [Fe/H] relation on the dif- 
ferent temperature scales produces slope values that are both 
always significant at the level of 3cr and, strikingly, consistent 
with each other within Icr. The slope values are also the highest 
reported to date. The hypothesis that the "slope" in the A(Li) - 
[Fe/H] relation might be due to the specific Teff scale chosen can 
thus be safely rejected. 

7.2. The meltdown of the Spite plateau: slope or scatter? 

Three different Ha-based Teff scales, as well as the totally inde- 
pendent IRFM scale, concur in indicating that the Spite plateau 
is disrupted below [Fe/H] ~ -3 (see Fig. 10). Close to that metal- 
licity, one observes a significant increase in the Li abundance 
scatter, which appears to act always towards lower abundances. 
In other words, while some rare stars persist at the plateau 
level even at very low metallicity (CS 22876-032 A, Gonzalez 
Hernandez et al., 2008), the vast majority exhibit some degree of 
Li "depletion" (with respect to the plateau value). 

Two stars in the sample exhibit anomalously low Li abun- 
dances. The star CS 22188-033 exhibits a mild Li depletion with 
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Fig. 10. Li abundance versus [Fe/H] for the four temperature estimates. 
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dances, while the red crosses represent ID LTE Li abundances with 
the NLTE corrections applied. The best-fit linear relation (as per Table 
8) is indicated by a gray line. A typical error bar of ±0.09 dex in Li 
abundance and the average [Fe/H] error bar are also displayed. 



a 3D NLTE A(Li) = 1 .66, while CS 22882-027 has no detectable 
Li doublet (see sect. 7.5). 

One of the much-debated results concerning the behavior of 
Li abundances in metal-poor halo dwarfs has been the reported 
existence of a correlation between A(Li) and [Fe/H], since it 
was first reported by Ryan et al. (1999). We investigated this 
by means of two different statistical tests. The Kendall's r rank- 
correlation test attempts to detect a (positive or negative) correla- 
tion, and has the fundamental strength of being non-parametric. 
In other words, it does not attempt to look for a specific rela- 
tion to fit the data. As seen above, Kendall's r-test quite strongly 
supports the existence of a correlation. 

The other obvious strategy we adopt is to fit the data with 
a linear function and see whether the slope found is statistically 
significant. This significance might be weakened if the data are 
indeed correlated, but the underlying relation is not linear In our 
case, again, the slope of the linear relation is significant at 3cr for 
all the considered temperature scales. This a finding, however, 
does not imply that the underlying "physical" relation between 
[Fe/H] and A(Li) is linear, as it would be, for example, if there 
was a constant Li production with increasing [Fe/H] . 
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Table 7. Parameters of the models in the 3D CO^BOLD and ID LHD grids 
used in the 3D NLTE Li abundances, and in the computation of NLTE 
corrections. 



Teff 


logg 


[Fe/H] 


ACLi)" 


EW 


K 


cgs 








pm 


5472. 


4.00 


-2.0 


0.90- 


2.10 


5.71 


- 69.71 


5479. 


4.50 


-2.0 


0.90- 


2.10 


5.94 


- 72.80 


5505. 


3.50 


-2.0 


0.90- 


2.10 


4.98 


-61.94 


5846. 


4.00 


-3.0 


1.30- 


2.50 


6.67 


- 76.04 


5856. 


4.00 


-2.0 


1.30- 


2.50 


6.79 


- 77.83 


5861. 


3.50 


-2.0 


1.30- 


2.50 


6.28 


- 72.67 


5923. 


4.50 


-2.0 


1.30- 


2.50 


6.17 


-73.16 


5924. 


4.50 


-3.0 


1.30- 


2.50 


5.76 


- 68.94 


6269. 


4.00 


-3.0 


1.30- 


2.50 


3.29 


- 43.33 


6272. 


4.50 


-3.0 


1.30- 


2.50 


3.32 


- 43.83 


6278. 


4.00 


-2.0 


1.70- 


2.90 


8.05 


- 86.35 


6287. 


3.50 


-2.0 


1.70- 


2.90 


7.42 


- 81.33 


6323. 


4.50 


-2.0 


1.70- 


2.90 


7.71 


- 84.35 


6408. 


4.00 


-3.0 


1.30- 


2.50 


2.68 


- 36.35 


6533. 


4.50 


-2.0 


2.10- 


3.30 


13.55 


- 119.59 


6534. 


4.00 


-2.0 


2.10 - 


3.30 


12.87 


- 118.58 


6556. 


4.50 


-3.0 


1.30 - 


2.50 


2.16 


- 30.20 



Mimmum and maximum A(Li) covered in the COG 
* Minimum and maximum EW for the Li doublet in the COG 

Table 8. Kendall rank correlation probability, intercepts, and slopes of 
the linear fit and standard deviation of the slope for the plots shown in 
Figs. 12 and 10. The star CS 22188-0033 has been excluded from the 
fit. 



Parameter 

set 



correlation 
probability 



linear fit 
intercept 



linear fit 
slope 



A(Li) vs. TeB 



BA 
ALl 
IRFM 
3D 



0.978 
0.962 

0.993 
0.990 



-3.038± 1.542 
-7.917±10.83 

-2.133±1.448 
-2.910±1.898 



A(Li) vs. [Fe/H] 



BA 
ALl 
IRFM 
3D 



1.000 
0.999 
0.998 
0.999 



3.099±0.245 
2.942±0.245 
3.047±0.268 
2.948±0.248 



0.335±0.080 
0.258±0.085 
0.291 ±0.092 
0.274+0.083 



To shed more light on the issue, in Fig. 11 we plot the residu- 
als of the best-fit A(Li) versus [Fe/H] relation, as listed in Table 

8. An increase in the scatter below [Fe/H] 2.8 was already 

visually apparent in Fig. 10, and remains clearly recognizable 
in Fig. 1 1 once the best-fit linear relation is subtracted. To pro- 
vide quantitative estimates of the level of scatter, we divided the 
sample into two in terms of metallicity, a metal-richer subsample 
including the 8 stars with [Fe/HJso > -2.8, and a metal-poorer 
sub-sample including the 19 stars below that threshold. The star 
CS 22188-033 is plotted in the figure, but it has not been con- 
sidered in this computation. We then computed the dispersion 
in the residuals of the two subsamples for the four temperature 
scales: o"hi,BA = 0.04 dex, cti^ba = 0.10 dex; ctiu^ali = 0.05 
dex, crio,ALi = 0.08 dex; ctm.irfm = 0.02 dex, ctio.irfm = 0.10 
dex; cTj,; 3o - 0.05 dex, ctio^d = 0.09 dex. For every temperature 
scale, the scatter in the residuals is about twice as large (or more) 
below [Fe/H]=-2.8 than above. It is thus clear that the tight, flat 
relation that is known as the Spite plateau develops both a tilt 
and a significant scatter at low metallicities. Once again, it is 
remarkable how the level of scatter appears independent of the 
assumed temperature scale. 
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Fig. 11. Residuals of the best fit of A(Li) vs. [Fe/H] listed in Table 8. 



We note that the true Li doublet EW does not change much 
with metalUcity, since in general A(Li) does not vary by more 
-than 0.2 dex. The quality of the Li doublet measurement is thus 
roughly constant across the whole metallicity range. The in- 
crease in scatter thus cannot be attributed to the dechning quality 
of the measurements. On the other hand, Fei and Fen Unes do 
become weaker with metallicity, which lowers the quality of the 
gravity estimation. Since the Ha line is quite gravity sensitive, 
inaccurate gravities reflect directly on Ha-based Tetf estimations, 
and thus on A(Li). On the other hand, the IRFJVI temperature 
scale is totally insensitive to this effect, and yet shows the largest 
increase in the scatter of its residuals, and a low-metallicity scat- 
ter equal to those of the Ha-based Tetr scales. This reinforces our 
impression that the increase in the A(Li) residual scatter should 
indeed be real. 



7.3. Plateau placement 

As a consequence of what is said above, it hardly makes sense 
to provide an average value for A(Li) in our stars. One might 
still try, however, to determine the position of the plateau for the 
more metal-rich stars of the sample, which still appear to fall 
onto it. Every operation of this kind is somewhat arbitrary, since 
there is no clear-cut transition between the plateau at higher Fe 
content and the sloping / dispersed distribution at low metallicity. 
We thus decided to employ the 9 stars whose metalhcity is equal 
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Fig. 12. Li abundance versus effective temperature for the four temper- 
ature estimates. Symbols are the same as in Fig. 10. Typical error bars 
of ±0.09 dex in the Li abundance and ±130 K in Teff are also displayed. 
The best-fit linear relation as per Table 8 is indicated by a gray line. 



or greater than -2.8 in the 3D scale, and consider their A(Li) 
and dispersion as being representative of the Spite plateau. The 
resulting values are <A(Li)) = 2.164 + 0.059 in the BA scale, 
(A(Li)) = 2.261 + 0.053 in the ALI scale, <A(Li)) = 2.264 + 
0.044 in the IRFM scale, and (A(Li)> = 2.199 + 0.086 using 
A(Li)3D NLTE together with the 3D temperature scale. 

7.4. The A(Li) - Teff correlation 

In Fig. 12, as well as in Table 8, a positive slope of Li abun- 
dance with effective temperature is present in all the temperature 
scales. It is, however, remarkable that this slope is driven by the 
coolest stars. Adopting our 3D temperature scale, if we select 
only the stars hotter than 6250 K (22 stars), the slope essentially 
vanishes, the probability indicated by Kendall's t drops to 89%, 
and a parametric test does not detect any slope. 

Given the very small population of this temperature range 
in our sample, it is possible that we just missed any undepleted 
cool object. On the other hand, this might suggest that the de- 
cline (caused by convection) usually seen for stars cooler than 
~ 5700 K at higher metallicities may set in at higher tempera- 
tures for EMP stars (although observationally the opposite seems 
to be true, see Boesgaard et al., 2005). 




517.22 517.24 517.26 517.28 517.30 

wavelength (nm) 

Fig. 13. The Mgi 517.268 nm line for HE 1148-0037, in spectra taken 
at JD 2453787.6803 (black continuous line) and JD 2453823.5443 (red 
dashed line), separated by 35.86 days. The spectra have been normal- 
ized, and Doppler shifted to bring the primary component to a rest- 
frame wavelength. The binarity is readily visible, as well as the varia- 
tion in the components' separation. 



We emphasize, however, that the depletion of the cool end of 
the sample is not driving the A(Li)-[Fe/H] correlation: removing 
the aforementioned five cool stars has a negligible effect on this 
result. On the 3D scale, the Kendall's t correlation probability 
of A(Li)3D,NLTE with [Fc/H] passes from 0.999 to 0.998, while 
the slope of the linear fit goes from 0.274±0.083 to 0.253+0.086 
when these 5 stars are removed. This is because the cool stars, 
while appearing to be all Li depleted, are evenly distributed in 
metallicity between [Fe/H] — 3.2 and -2.8. On the other hand, 
removing these stars would affect the detected increase in the 
scatter at lower metallicities: applying the same residual analysis 
mentioned in Sect. 7.2, but now on the hot stars only, would 
yield crj,; 3^ ^ot = 0.04 dex (was 0.05 with the full sample), and 
CTio.sD.hot = 0.06 dex (was 0.09 with the full sample). 

7.5. CS 22882-027 and HE 1 148-0037 

Two additional stars were included in the original sample, but 
Li abundance has not been computed for them, for different rea- 
sons. We elaborate briefly on these objects. 

Atmospheric parameters and metallicity were computed for 
CS 22882-027 (see Table 4), but the star has no detectable Li 
doublet (see Fig. 1). Given a S/N~80, as measured in the Li 
doublet range, a typical line FWHM of 0.033 nm, and a pixel 
size of 0.0027 nm, the Cayrel (1988) formula predicts that a 
Li doublet of EW=0.563 pm would be measured at 3cr confi- 
dence (lcr=0.187 pm). Employing our fuctional fit, this leads to 
an upper limit of A(Li)3D,NLTE ^ 1.82, assuming an EW=0.563 
pm and A(Li) 3d,nlte <1.34, assuming EW=0.187 pm. We have 
a single-epoch spectrum for this star that shows no sign of a 
double-line system. 

The star HE 1 148-0037 was also originally included in the 
sample, but immediately set aside, since from visual inspection 
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Table 9. Barycentric radial velocities for the two components of the 
binary system HEl 148-0037 as measured from the spectra available to 



date 
UT 



hour 
UT 



JD 

UT 



^ rad ^ rad 

km/s km/s 
primary secondary 



8 May 2003" 
27 Feb. 2005* 
21 Feb. 2006 
29 Mar. 2006 



11:43:06 
10:33:36 
04:19:40 
01:03:46 



2452767.9882 
2453428.94 
2453787.6803 
2453823.5443 



-36.16±1 -3.35±1 

-10.88±0.16 
-9.15±1 -22.37±1 
-6.90±1 -22.56±1 



a HERES spectram 

h Aoki et al. (2009) single radial velocity 



of the two available spectra it turned out to be a double-lined 
binary system. As shown in Fig. 13, the two spectra, separated 
by 35.86 days, clearly exhibit evidence of the double-Une sys- 
tem, as well as readily recognizable variation in the separation 
between the two line systems. We were able to retrieve 3 spec- 
tra of HE 1 148-0037, and measure radial velocities for the two 
components, which are given in Table 9, by means of cross- 
correlation against a synthetic template (reff =6000 K, log g=4.0, 
[Fe/H]=-2.0). In our two spectra (2006-02-21 and 2006-03-29), 
cross-correlation was computed in the 490nm-570nm range. We 
also had the lower resolution, blue-range only HERES spectrum, 
which was used to derive the Vrad after masking all the broad 
hydrogen lines. Internal errors of the radial velocity estimate 
were evaluated by performing a Monte Carlo test on a sample 
of 50 simulated binary star spectra with noise added to emulate 
a S/N=90, the component separation and resolution being equiv- 
alent to that of the March 29, 2006 observation. The test inferred 
an average error of 0.072 km/s for the primary component and 
0.246 km/s for the secondary. These values are representative 
of the internal errors in the cross-correlation procedure, but are 
surely dominated by the spectrograph zero-point cahbration un- 
certainty, which was not taken into account particularly well, be- 
cause very precise radial velocities were not among the goals of 
this study. We thus list in Table 9 an estimated total error of 1 
km/s for all our observations, which we take to be representative 
of the overall systematics of these measurements. 

On the other hand, Aoki et al. (2009) report AN uncertainty 
given by the internal scatter in the radial velocities obtained 
from different Unes; they do not take into account the system- 
atic uncertainties in the wavelength calibration, hence the much 
lower value of uncertainty. In Tab. 9, we kept the value they pro- 
vide, but we propose that the true uncertainty is again close to 
1 km/s. If we consider their measure to be representative of the 
primary component radial velocity, it reproduces well the three 
other measurements we present. 

7.6. Comparison with ottier results 

Our most significant overlap is of course with the Bonifacio et 
al. (2007) sample, of which the present study represents a con- 
tinuation. Of the 19 stars in Bonifacio et al. (2007), 17 were 
reanalyzed here (the two remaining stars, BS 16076-006 and 
BS 16968-061 turned out to be subgiants and have thus been 
dropped). Bonifacio et al. (2007) determined Teff by fitting the 
Hawings fits with profiles synthesized using the Barklem et al. 
(2000a,b) self-broadening theory. Their temperature scale thus 
closely resembles our BA scale, but the Ha gravity sensitiv- 
ity was not taken into account in that study - log g=4.0 was 
assumed in computing the profiles. The effect is clearly visi- 
ble in Fig. 14. The upper panel of this figure plots the tem- 



perature difference, Teff (this work)-reff(Bonifacio et al. 2007), 
against the Bonifacio et al. (2007) gravity estimate. It is im- 
mediately evident that, for stars around log g=4.0, the tempera- 
4ure difference approaches zero. For stars at higher gravities, our 
Teff estimate is below the Bonifacio et al. (2007) value by up to 
200 K. This reflects the behavior presented in Fig. 3: when a pro- 
file is fitted with a synthetic grid computed for a gravity that is 
underestimated, this leads to an overestimated T^g. We note that 
Fig. 14 does not tell the entire story. As we do here, Bonifacio 
et al. (2007) estimated log gby enforcing Fe ionization equilib- 
rium, so the derived gravity values are not identical in this work 
and Bonifacio et al. (2007). High-gravity and low-gravity stars 
do, however, retain their approximate placement in both cases, 
although the gravity span can be somewhat stretched by the Teff 
bias. The effect of the T^g difference on [Fe/H] is shown in the 
middle panel of Fig. 14, while in the lower panel the difference 
between A(Li) for the same stars is plotted, considering here the 
LTE values (to eliminate the effect of the marginally different 
NLTE corrections applied in the two works). The results shown 
are to be expected, given the strong T^g sensitivity of the Lii 
670.8 nm doublet: the current A(Li) is higher by up to about 0.1 
dex for low-gravity stars, while it is lower by roughly the same 
amount for high-gravity stars. On the other hand, it is easy to 
see how the discrepancy will only marginally affect a linear fit 
of A(Li) versus [Fe/H]: the stars are displaced roughly along a 
1:1 diagonal in the A(Li) - [Fe/H] plane. 

The star LP 815-43 is the only object that overlaps with the 
Asplund et al. (2006) sample. In that work, effective tempera- 
ture is measured again by fitting Hor wings with a set of syn- 
thetic profiles. The details of the fitting procedure differ some- 
what and the synthetic profiles are computed from MARCS models 
using the BSYN synthesis code. The Barklem et al. (2000a,b) self- 
broadening theory is assumed here for Ha, so again the BA scale 
is the one to be used in the comparison. The Hor gravity sensi- 
tivity is taken into account here, and the derived T^g is quite 
close to our value (6400 K vs. 6453 K in this work). Asplund 
et al. (2006) determine metallicity from Fe ii lines, and gravity 
from Hipparcos parallaxes, but again the values do not differ 
much from our results (log g=4.17, Vturb=1.5, [Fe/H]=-2.74). 
The residual 0.14 dex offset in metallicity is in good agreement 
with the 0.2 dex offset detected by Bonifacio et al. (2007) be- 
tween their metallicity scale and that of Asplund et al. (2006). 
The hthium abundance is again quite close to our result: from 
the 670.8 nm doublet, they derive A(^Li)=2.16, while our value 
is 2.23. The difference is fully accounted for once the T^g effect 
is considered (0.03 dex) as well as the already known, albeit still 
unexplained, 0.04 dex bias between A(Li) as derived by means 
of BSYN and turbospectrum (Bonifacio et al., 2007). 

The same star is also in common with Hosford et al. (2009). 
That work uses excitation equilibrium to estimate T^^g, a method 
that is not directly comparable with any of our temperature scale. 
The authors derive two parameter sets, one assuming the star 
belongs on the main sequence, and another assuming it is a 
subgiant. In the first case (reff=6529, log g=4.40, Vturb=l-4, 
[Fe/H]=-2.61), they obtain a temperature that is only about 
50 K cooler than for our 3D and ALI scale, but a higher grav- 
ity and metalUcity. In the second case (reff=6400, log g=3.80, 
Vturb=l-4, [Fe/H] =-2.68), they derive a temperature 50 K cooler 
than for our BA scale, the same gravity, but again to a metallic- 
ity 0.2 dex higher. Their measured Li doublet equivalent width 
is about 0.2 pm smaUer than that we find, and in both cases they 
derive an A(Li) that is about 0.1 dex smaUer than ours. This is 
consistent with the expected combined effect of the difference in 
both Tgg and the Li doublet EW. 
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Fig. 14. Comparison between the results of this work and of Bonifacio 
et al. (2007) for the 17 stars in common. In the upper panel, the differ- 
ence between our and the Bonifacio et al. (2007) T^g determination 
is plotted against the value of log g in Bonifacio et al. (2007). In the 
center panel, we show [Fe/H] difference, and in the lower panel, the 
difference in the A(Li) LTE. For our results, BA temperature scale is 
used. 



Three stars are in common with the Aoki et al. (2009) sam- 
ple, but we determined A(Li) for only two of them, CS 22948- 
093 and CS 22965-054. Aoki et al. (2009) employ Hor- as well 
as I^-wing fitting to determine T^^g on the basis of MARCS mod- 
els and using the Barklem et al. (2000a,b) self-broadening for 
Ha. Once more, our BA temperature scale is the one most ap- 
propriate for a comparison. The Ha gravity sensitivity is taken 
into account. Surface gravities are estimated by comparing with 
isochrones as well as by evaluating the Fe i-Fe ii ionization equi- 
librium. The parameter values they derive for CS 22948-093 
are in close agreement with those we derive (reff(Ha)=6320 K, 
reff=6380, log g=4.4, Vturb=1.5), while [Fe/H] is 0.13 dex lower 
at -3.43. The derived value of A(Li)LTE=l-96 is inexcellent 
agreement with our value of 1 .935. CS 22965-054 shows a more 
significant discrepancy in T^s. Their value of reff(Ha)=6390K is 
remarkably higher than our, while Tgff(H/3) is much closer to our 
value. Since the adopted temperature is the average of the two, 
the final Teff is ultimately just 56 K hotter than our value. Their 
derived surface gravity is also very close in value (log g=3.9), 
while Vturb is the same. However, Aoki et al. (2009) derive the 



same lithium abundance as we do, A(Li)LTE=2.16, due to their 
lower measured value of EW (2.03 pm) for the Li doublet. 

The third star in conmion with Aoki et al. (2009) is HE 
1 148-0037, for which Aoki et al. (2009) do not appear to have 
noticed its binarity. On close inspection, the single spectrum 
they employed (JD 2453428.94) shows signs of line asymme- 
try (Aoki 2009, priv. comm.), but it appears to have been taken 
quite close to conjunction. In contrast, as seen in Sect. 7.5, the 
three spectra we have of this star all show quite clearly the two 
line systems. 

We observed No star in common with Melendez & Ramirez 
(2004). A comparison with their sample is nevertheless inter- 
esting due to the extension of their sample to very low metal- 
licities, despite with a limited number of stars (a total of 10 
stars were analyzed below [Fe/H]=-2.5, 4 at or below [Fe/H]=- 
3). Melendez & Ramirez (2004) detected no slope in the Spite 
plateau, for which they advocated a high value of (A(Li))=2.37. 
As noticed in Gonzalez Hernandez & Bonifacio (2009), the T^^g 
scale adopted by Melendez & Ramirez (2004) is systematically 
hotter than the one we employ for metal poor dwarfs, on average 
by 87 K. This explains half of the descrepancy between our av- 
erage IRFM plateau placement and their own, and can account 
in principle for their failure to detect the slope, assuming their 
temperature scale and ours diverge progressively at low metal- 
licities. 

Lithium abundances for two extremely metal poor stars (HE 
0233-0343 and HE 0945-1435) were recently presented by 
Garcia Perez et al. (2008). Both stars show extremely low Fe 
content ([Fe/H] — 4), but probably because of the weakness of 
Fen lines, the estimation of gravity is uncertain. This affects 
the determinations of both the evolutionary status (either MS 
or early SGB) and T^s, which is derived from Ha wing fitting 
in a way similar to that used with our BA scale. The stars ap- 
pear to be fairly cool, 6000K<reff<6250K, which would place 
them among the "cool stars" of our sample as described in Sect. 
7.4, and both objects show significantly depleted Li, A(Li)~1.8. 
Owing to the uncertainty of the parameters determination we did 
not include these stars in Fig. 15. 

8. Possible biases 

8.1. Binary stars 

Two potential biases can, in principle, be responsible for produc- 
ing systematically low Li abundances, and a trend of A(Li) with 
metallicity. The first one is of course the presence of undetected 
binaries, for which the veiling by the secondary star will sys- 
tematically reduce the EW of the lines of the primary, leading 
to an underestimate of both metallicity and Li abundance. The 
true impact of this effect is difficult to evaluate, mainly because 
little is known about the fraction, and mass-ratio distribution, of 
binaries at low metallicities. Duquennoy et al. (1991) report, for 
G dwarfs in the solar vicinity, a fraction of 44% of stars hav- 
ing a companion with q = M2/M1 > 0.1, about 1/3 of which 
have q > 0.5. Latham et al. (2002, and references therein) found 
that the halo binary population does not differ significantly from 
the disk one, although we note that they did not explore signifi- 
cant numbers of stars with metallicities as low as the stars in our 
present sample. 

It nevertheless seems unlikely that undetected binaries pol- 
lute our sample significantly. We checked for binarity by inspect- 
ing the Mg b triplet lines (e.g., see Fig. 13). Our spectra have a 
typical S/N~100 or higher, and Mg b lines have a typical EW 
of 10 pm. A line with central residual intensity of 0.95 would 
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be detected at least at the 5cr level in this typical spectrum, and 
have a typical EW of 0.8 pm. As per Gonzalez Hernandez et al. 
(2008), to reduce a 10 pm line to 0.8 pm, a ratio of the continua 
fluxes of about 11.5 is needed. If we roughly assume that the 
total luminosity scales accordingly, this corresponds roughly to 
q = 0.5 (since L oc M^-^ on the main sequence, see Kippenhahn 
& Weigert 1990). A similar flux ratio in the Li doublet range 
would lead to a correction of the Li doublet EW for the primary 
star of about 8%, corresponding to 0.03 dex in A(Li). In other 
words, every binary star requiring significant veiling correction 
on the primary spectrum would also be promptly detectable be- 
cause of the double Une system. This system could only stay 
undetected if the radial velocity separation of the two stars was 
quite small at the moment of the observation(s), so that the two 
fine systems remained blended. 

In addition to the above, if we take the figures of Duquennoy 
et al. (1991) at face value, find that about 13.5% of the bina- 
ries are characterized by a significant veiling of the primary (i.e. 
q > 0.5). Our original sample comprised 30 stars, which implies 
that there are 4 expected "significant binaries". Three stars have 
already been rejected from the sample, one of them (HE 1148- 
0037) being indeed a binary. Two more (CS 22882-027 and CS 
22188-033) exhibit significant lithium depletion or no Li dou- 
blet at all, and were excluded from all statistical analyses. They 
are clearly the most Ukely candidates to be binaries, albeit nei- 
ther one shows a double line system^. One could thus expect one 
more "disguised" binary to be biasing the sample. While this is 
quite possible, it would hardly influence any of our results. 

8.2. 3D NLTE effects on Fe ionization equilibrium 

The second problem relates to the use of Fei-Fen ionization 
equiUbrium to estimate gravity. From preliminary computations, 
it appears that 3D corrections of Fe fines with excitation poten- 
tials of the same order as employed in the present work could be 
quite large at low metallicities for stars similar to those that we 
study. Moreover, corrections for Fe i appear to be negative (of 
about 0.2 dex), while they are positive (about 0. 1 dex) for Fe n 
fines, for a reff=6500, log g=4.5, [Fe/HJ=-3.0 star. The phe- 
nomenon is mainly caused by the overcooling that 3D treatment 
produces in the outer layers of atmospheres at low metalUcities, 
and appears to be of similar magnitude at [Fe/H]=-2 (due to 
the stronger saturation of Fe lines, which drives their contribu- 
tion function to higher layers), but would most likely disappear 
above. If taken at face value, a 0.3 dex Fe i-Fe n imbalance would 
lead to an overestimate of log g of about 0.5 dex when analyzed 
using ID LTE models (as is our case regarding metallicity and 
gravity estimation). We do indeed find higher gravities than ex- 
pected from evolutionary tracks. On the other hand, we do not 
currently have a 3D NLTE spectrosynthesis code for iron; we 
are thus unable to account for NLTE effects, which are likely to 
counterbalance the 3D effect because of over-ionization occur- 
ring in the upper layers where overcooling is present in 3D mod- 
els. A similar mechanism is indeed active for Li, whose 3D LTE 
abundance derived by CO^BOLD-LINFORBD is about 0.2 dex be- 
low the corresponding ID NLTE value, while the 3D NLTE one 
is essentially indistinguishable from the ID NLTE result (Figs. 
10 and 12). If some degree of imbalance remains after NLTE 
is taken into account, and is metalUcity sensitive, this might in- 



' We have another UVES spectrum of CS 22188-033, taken at a dif- 
ferent epoch, which does not show radial velocity variations with re- 
spect to the one used in the present work, nor signs of a double line 
system. 



troduce a bias when T^s is determined by Ha-wing fitting. It 
is indeed intriguing to note how lower gravity estimates lead to 
higher temperatures, and as a consequence, (somewhat) higher 
the Li abundances. On the other hand, the IRFM temperature 
scale should be immune to this problem, A(Li) being quite in- 
sensitive to log g itself, and the IRFM-based analysis inferring 
abundances similar to the Ha-based estimate. 



9. Conclusions 

We have presented the largest sample to date of Li abundances 
for EMP halo dwarf stars (27 abundances and one upper limit), 
including the largest sample to date below [Fe/H] = -3(10 abun- 
dances). Lithium abundance determination is highly sensitive to 
biases in the effective temperature scale, and we have tried to ac- 
count for this using four diff'erent temperature estimators. In an 
additional effort to accurately represent the stellar atmospheres 
of the sample stars, 3D, time-dependent, hydrodynamical atmo- 
sphere models have been used to determinine our preferred Hor- 
based temperature scale, and a detailed 3D NLTE spectrosynthe- 
sis has been applied to the determination of lithium abundance. 
Both these techniques have been employed here for the first time, 
to our knowledge, in the analysis of EMP stars. This has also 
allowed us to develop a useful fitting formula allowing one to 
derive A(Li)3D,NLTE directly as a function of EW, Teff, log g, and 
[Fe/H] for EMP turn-off and early subgiant stars (see Appendix 
B). 

The first obvious conclusion of this work is that we have con- 
firmed what was merely suggested by the analysis of Bonifacio 
et al. (2007), and previous works, that at the lowest metallicity 
there is sizable dispersion in the Li abundances and that there 
is a trend of decreasing Li abundance with decreasing metallic- 
ity. We have also shown that these two conclusions do not de- 
pend on the adopted temperature scale, as suggested by Molaro 
(2008). The results hold, qualitatively, using both IRFM temper- 
atures and Ha temperatures, regardless of the broadening theory 
adopted and irrespective of the use of either ID or 3D model 
atmospheres. Quantitatively, the results differ in the mean level 
of the Li abundance, while the slopes in the A(Li) versus [Fe/H] 
relations agree within errors. None of the temperature scales in- 
vestigated produces a "flat" Spite plateau over the full range in 
[Fe/H] (see Table 8). 

Our results are in substantial agreement with those of Aoki et 
al. (2009). While these authors do not detect a slope with either 
effective temperature or metallicity, this happens simply because 
of the smaU extent of their sample in both these parameters. On 
the other hand, they do point out that their sample has a lower Li 
abundance than that observed at higher metalUcities. 

The picture outlined by the aforementioned results acquires 
more significance, once we place it in a broader context among 
the latest studies regarding lithium in EMP stars. Figure 15 com- 
pares our results with those of three investigations employing 
compatible temperature scales. In this figure, open blue circles 
represent stars from the Asplund et al. (2006) sample, red tri- 
angles from the Aoki et al. (2009) data, and the two magenta 
squares the two components of the double-lined binary sys- 
tem CS 22876-032 (Gonzalez Hernandez et al., 2008, flie fiUed 
square corresponds to the primary star). Our data are repre- 
sented as black diamonds (the results of the BA scale are shown, 
for compatibility with the temperature scales used in the other 



L. Sbordone et al.: The metal-poor end of the Spite plateau 



17 



< 



2.6 










o ; 


2.4 












2.2 






o 


A R 


^ O O CSCr 
- 


2.0 




^/^O o 

o 








1.8 




□ o 








1.6 




o 

1 .... 1 


1 




1 



-3.5 -3.0 -2.5 -2.0 -1.5 

[Fe/H] 



Fig. 15. A unified view of A(Li) vs. [Fe/H] from some studies for which 
a common temperature scale can be assumed. Blue circles, Asplund et 
al. (2006) data, red triangles, Aoki et al. (2009) data, magenta squares, 
CS 22876-032 from Gonzalez Hernandez et al. (2008), filled symbol 
primary star, open symbol secondary star. Black diamonds, this work, 
BA temperature scale. Dot-dashed gray line, best linear fit to Asplund 
et al. (2006) data, continuous dark gray line, best fit to our data. Typical 
error bars for our data are displayed. 



three works) . The best linear fit to our data is shown as a dark 
gray solid line, while the best fit to Asplund et al. (2006) data 
(A(Li)=2.409-H0.103[Fe/H])is shown by a dot-dashed gray line. 
The Asplund et al. (2006) Li abundances are increased here by 
0.04 dex to account for the known offset already mentioned in 
Sect. 7.6, and their metallicty is decreased by 0.2 dex to corre- 
spond to the metallicity-scale offset detected by Bonifacio et al. 
(2007). It is now even more evident that the Spite plateau does 
not exist anymore at the lowest metallicity, and is replaced by an 
increased spread of abundances, apparently covering a roughly 
triangular region ending quite sharply at the plateau level. This 
region appears here to be populated in a remarkably even man- 
ner; at any probed metallicity some star remains at, or very close 
to, the Spite plateau level, but many do not. The rather different 
slopes of the best-fit relations in Asplund et al. (2006) and in 
this work appear to be the obvious consequence of fitting two 
subsamples covering different metallicity regimes. This could 
provide also an explanation for the numerous claims, starting 
from Ryan et al. (1999), of a thin, but tilted Spite plateau. From 
this view, the difference was produced simply because the tail of 
these samples had been falling in the low-metallicity "overde- 
pletion zone" as we have been able to discern more clearly. 

We are not aware of any theoretical explanation of this be- 
havior. After the measurements of the fluctuations of the CMB 
made it clear that there is a "cosmological lithium problem", i.e., 
the Li predicted by SEEN and the measured baryonic density is 
too high with respect to the Spite plateau (by about 0.6 dex for 
our sample), there have been many theoretical attempts to pro- 
vide Li-depletion mechanisms that would reduce the primordial 
Li to the Spite plateau value in a uniform way. Our observations 
now place anadditional constraint on these models - below a 
metallicity of about [Fe/H] = -2.5, they should cause a disper- 
sion in Li abundances and an overall lowering of A(Li). 

If Li depletion from the WMAP-prescribed level were 
to happen in the stellar envelopes of very metal-poor stars. 



* Gonzalez Hernandez et al. (2008) derived Tcff from photometry and 
isochrones, but a cross-check with Ha profiles computed in ID with 
Barklem et al. (2000a) broadening confirmed the result. 



the mechanism would have to be remarkably metallicity in- 
sensitive to account for the thin, flat plateau observed be- 
tween [Fe/H]=-2.5 and -1. And yet, the same phenomenon 
must become sharply metallicity sensitive around and below 
[Fe/H] =-2.5, i.e., precisely where metallicity effects on the at- 
mospheric structure are expected to become vanishing small. 

We are tempted to imagine that two different mechanisms 
may need to be invoked to explain the production of the Spite 
plateau for stars with [Fe/H] > -2.5, and of the low-metallicity 
dispersion for stars with [Fe/H] < -2.5. One could envision such 
a two-step process as follows: 

1 . Metal-poor halo stars are always formed at the Spite plateau 
level, regardless of their metallicity. Whether the plateau rep- 
resents the cosmological Li abundance or is the result of 
some primordial uniform depletion taking place before the 
star formation phase is immaterial in this context. 

2. A second phenomenon, possibly related to atmospheric dif- 
fusion, becomes active around [Fe/H]=-2.5 and below, de- 
pleting Li further in the atmosphere of EMP stars. This phe- 
nomenon, aside from the metallicity sensitivity, would ex- 
hibit different star-to-star efficiency, being possibly depen- 
dent on additional parameters, such as stellar rotation or Teff . 
Its efficiency must in any case be higher for more metal-poor 
stars. 

In this scenario, the "primordial" plateau would be preserved 
above [Fe/H]~ -2.5, but below that metallicity, a systematic 
"leakage" of stars towards lower A(Li) would take place, more 
effectively for more metal-poor stars, but naturally scattered due 
to the sensitivity to parameters other than [Fe/H]. This scheme 
would have a number of advantages. First of all, it would natu- 
rally explain our observations, "mimicking" a slope in A(Li) ver- 
sus [Fe/H], but with increased scatter at low [Fe/H]. It would also 
explain why, while the scatter in A(Li) increases at low metallic- 
ities, not a single star in this metallicity regime has been found 
to lie above the Spite plateau level. It would then be consistent 
with a small number of stars remaining close to the plateau at 
any metallicity (e.g., CS 22876-032 A, Gonzalez Hernandez et 
al., 2008, filled magenta square in Fig. 15); in these objects, 
the depletion process would be somehow inhibited. Finally, at- 
tributing the extra depletion to atmospheric diffusion / settling 
would not require a physical "conspiracy" capable of producing 
exactly the same depletion level regardless of metallicity, stel- 
lar rotation, gravity, or effective temperature, as is often invoked 
when diffusion is used to explain the Spite plateau. 

The nature of what we refer to above as the "second phe- 
nomenon", the one responsible for the departures from the Spite 
plateau below [Fe/H] - -2.5, is perhaps the most intriguing. 
Above, we have proposed some kind of photospheric settling 
mechanism, but one could as well envision a chemical evolution 
scenario, on the basis of some gas pre-processing with Li deple- 
tion (fl la Piau et al. 2006) - while it may not be able to account 
for the entire WMAP-Spite plateau discrepancy, this mechanism 
could easily account for the mild (0.2-0.4 dex) departure from 
the plateau observed at lower metallicities. Moreover, this mech- 
anism would naturally produce a spread of abundances as a con- 
sequence of the local level of gas pre-processing. 

There are hints that the recently discovered ultra-faint dwarf 
galaxies (uFdg) might have been the source of the bulk of the 
EMP stars now found in the halo of the Milky Way (Tolstoy et 
al., 2009, and references therein). If this were indeed the case, a 
sizeable fraction of our sample could have formed in uFdg sys- 
tems, possibly more so for the most metal-poor objects. It has 
been suggested (Komiya et al., 2009) that the paucity of stars 
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below [Fe/H] = -3.5 may be due to the onset of self-pollution 
in the primordial mini-halos when they started to merge to form 
larger structures, and ultimately the halo. One could then envi- 
sion that gas reprocessing could have fairly significantly altered 
the Li abundances in the heavily dark-matter dominated cores of 
these sub-halos, but become progressively negligible when they 
merged, and their ISM mixed more and more completely with 
gas of pristine Li abundance. Sub-plateau stars may then origi- 
nate from star formation that occurred in the sub-halo cores, or 
while mixing with pristine gas that progressively diluted the ef- 
fect of the reprocessing. 

In any case, the main drawback of this two-phenomena sce- 
nario is to leave the WMAP-Spite plateau discrepancy unex- 
plained. It also implies that any depletion from WMAP-based 
primordial A(Li) should have taken place before the currently 
observed stars formed. Non-standard primordial nucleosynthe- 
sis thus remains entirely viable, until the contrary is proven. 

Looking retrospectively at the history of the abundances of 
light elements and SBBN, one concludes that if the Spite plateau 
is the result of significant uniform Li depletion among metal- 
poor stars, one faces a formidable case of cosmic conspiracy. It 
is necessary to admit that the Li has been depleted quite precisely 
to a level consistent with SBBN production and also consistent 
with the abundances of the other measurable light elements. Had 
Li been depleted to the level of A(Li)=1.8, that is, below the 
minimum allowed by SBBN, the Spite plateau might not have 
been interpreted as being related to primordial Li, right from the 
beginning. 

A caveat is in order when looking at Fig. 15: the Asplund et 

al. (2006) sample is biased, having been purposedly selected not 
to include objects significantly deviating from the Spite plateau. 
In fact, stars with varying degrees of Li depletion have long been 
known to exist (e.g., Charbonnel & Primas, 2005). On the other 
hand, warm (i.e., Teg- > 6000 K) dwarf stars deviating from the 
Spite Plateau appear to be quite rare at [Fe/H]>-3. Again, our 
sample is limited towards higher metallicities, which prevents us 
from properly quantifying this statement. A re-evaluation of the 
Li abundances above [Fe/H] =-2.5 is in order, using an analysis 
technique homogeneous with the present one, and applied to a 
sample for which a well-known selection technique has been ap- 
plied. It is plain to see that, had previous authors encountered Li 
abundance distributions as the one observed by us and Aoki et al. 
(2009), the concept of a lithium abundance plateau would have 
been short lived. Instead, there is a clear perception that mildly 
Li depleted stars do exist among MP stars, but are a rare occur- 
rence. We cannot easily identify any reason why all the previous 
studies should have been biased in a different way with respect 
to the present one, but until a coherent work, based on a well de- 
fined selection criterion, is extended to higher metallicities, all 
we are left with is just this perception. 

Two other Unes of future investigation are suggested by the 
present results: (i) increase the sample of hot (7'efr> 6250 K)) 
BMP stars ([Fe/H]< -3.0), and (//) increase the number of 
cool (reff< 6250 K) stars over the entire metallicity range. The 
present sample has seven hot BMP stars, five of which have 
A(Li) below the Spite plateau. Even if "standard" ZAMS Li de- 
pletion was setting in at higher than usual temperatures in BMP 
stars, these objects are unlikely to be affected by it. Statistics is 
still weak in this metallicity regime, and that the primary com- 
ponent of CS 22876-032 lies at the level of the Spite plateau, 
at a metalUcity of [Fe/H] = -3.6, clearly calls for a search for 
similar objects. At the same time, the cool end of our sample is 
presently too poorly populated to draw any definitive conclusion 
about whether these stars are experiencing convective depletion. 



which again, calls for an enlargement of the cool stars sample at 
low ([Fe/H]<-2.5) metallicities. 

To help understand the phenomena involved (since there 
could be several at work) that bring about the observed Li abun- 
dance pattern, another important issue is to assign an accurate 
evolutionary status to each star, that is, to confidently assess 
which stars are dwarfs, main-sequence tumofi", or subgiants. Our 
surface gravities are not sufficiently accurate for this purpose. 
The GAIA satellite will provide accurate parallaxes for all of the 
presently studied stars, thus allowing one to more clearly eluci- 
date this problem. 
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Appendix A: 3D non-LTE treatment of lithium 

The calculation of synthetic non-LTE Li I A 670.8 nm line pro- 
files from 3D hydrodynamical model atmospheres proceeds in 
two basic steps. In the first step, the code NLTE3D provides the 
departure coefficients biix,y,z,t) - ni{x,y,z,t)/n*(x,y,z,t), the 
ratio of non-LTE to LTE population number densities for each 
level i of the Li i model atom as a function of the geometrical 
position (x,y,z) in the 3D model atmosphere, and time (f) as 
sampled by a number of snapshots 20) selected to represent 
the characteristic temporal variation of the simulation. The de- 
parture coefficient of Lin is assumed to be 1, since lithium is 
essentially fully ionized in the stellar atmospheres of interest. 

In the second step, the departure coefficients are fed into the 
completely independent spectrum synthesis code LinforSD^, 
where they are used to compute the non-LTE line opacity and 
source function, and in turn the emergent intensity profiles as a 
function of (x, y, 6, (p, t), where the angles 6 and (f> specify the 
orientation of the line-of-sight in polar spherical coordinates. 
Finally, the emergent mean flux profile is obtained by horizon- 
tal, angular, and temporal averaging of the individual intensity 
profiles. 

For this paper, we used an 8 -level model atom of Lii to 
solve the statistical equilibrium equations, considering a total of 
1 1 bound-bound transitions. Details about the energy levels and 
line transitions are given in Tables A.l and A.2. The Einstein 
coefficients Aji provided by the MIST database are related to the 
Einstein coefficients Bij by 



C^ c2 gj 



(A.1) 



The corresponding oscillator strength fij is obtained from the 
relation 



gj MeC 
Jij - ^ji ~ 



niehcvij 



(A.2) 



The computationally most expensive part of solving the non- 
LTE problem is the calculation of the line-blanketed radiation 
field Jvix, y, z, t) at each grid point of the selected 3D models, 
which is needed to determine the photoionization rates for all 
atomic levels. This is done with a modified version of the ra- 
diation transport routines that are used in the CO^BOLD hydro- 
dynamical simulations for computing the radiative energy ex- 
change term VFrad(x,y,z). The solution of the radiative trans- 
fer equation is based on a Feautrier scheme applied to a set of 
long characteristics. The continuous opacities used in this con- 
text are computed with the routines lONDIS & OPALAM from 
the Kiel stellar atmosphere package'''. Line blanketing is taken 
into account by adding to the continuous opacity the opacity 



' http : //www . aip . de/~mst/Linfor3D/linfor_3D_raanual . 
pdf 

http : //www . aip . de/ ~ mst/Linf or 3D/linf or . pdf 
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Table A.2. Bound-bound transitions of the Li i model atom. Data are taken from the NIST database. 



transition 


lower 


upper 


transition 


A 


transition probability 


# 


level 


level 


configuration 


vacuum [A] 








1 


1 


2 


2s-2p 


6709.7 


3.72E-I-07 


8.49E-I-10 


7.532E-01 


2 


1 


4 


2s-3p 


3233.6 


1.17E-I-06 


2.99E-I-08 


5.502E-03 


3 


1 


7 


2s-4p 


2742.0 


1 .42E+06 


2.21E+08 


4.802E-03 


4 


2 


3 


2p-3s 


8128.6 


1.74E+07 


7.84E+09 


5.745E-02 


5 


2 


5 


2p-3d 


6105.3 


5.11E+07 


4.88E+10 


4.759E-01 


6 


2 


6 


2p-4s 


4973.1 


5.05E+06 


5.21E+08 


6.241E-03 


7 


2 


8 


2p-4d 


4604.1 


1.64E+07 


6.71E+09 


8.687E-02 


8 


3 


4 


3s-3p 


26887.1 


3.77E-I-06 


5.53E-I-11 


1.226E-I-00 


9 


3 


7 


3s-4p 


10795.1 


3.69E-I-03 


3.51E-I-07 


1.934E-04 


10 


4 


8 


3p-4d 


17550.0 


4.89E-I-06 


l.llE-i-11 


3.763E-01 


11 


5 


7 


3d-4p 


19281.0 


3.31E-I-05 


3.58E-I-09 


1.107E-02 



Table A.l. Energy levels of the Lii model atom. Data are taken from 
the NIST database. 



level 


configuration 


energy 


statistical 


# 




[Ryd] 


[eV] 


weight 


1 


2s 


0.0000000 


0.00000 


2 


2 


2p 


0.1358136 


1.84784 


6 


3 


3s 


0.2479204 


3.37313 


2 


4 


3p 


0.2818128 


3.83426 


6 


5 


3f 


0.2850726 


3.87861 


10 


6 


4s 


0.3190534 


4.34094 


2 


7 


4p 


0.3323350 


4.52165 


6 


8 


4d 


0.3337369 


4.54072 


10 



distribution functions (ODFs, 'big division', Vturb = 2 km/s) 
of Castelli & Kurucz (2003), including the Hi-H+ and Hi- 
ll i quasi-molecular absorption near 1400 and 1600 A, respec- 
tively. There is a slight inconsistency in the chemical compo- 
sition adopted for the calculation of the opacities: the continu- 
ous opacities are based on the solar abundances of Asplund et 
al. (2005), while the ODFs rely on the solar composition of 
Grevesse & Sauval (1998). To obtain the opacities for different 
metallicities, the solar abundances were scaled by a global fac- 
tor corresponding to the desired [M/H], with an enhancement 
of the a-elements by 0.4 dex below [M/H] = -0.5. A total of 
600 frequency points were used to obtain the Jy between A 925 
and 19 800 A. We checked that treating continuous scattering as 
true absorption does not introduce any significant changes in the 
resulting departure coefficients. 

Given the line-blanketed radiation field Jy{x, y, z, t), the pho- 
toionization rate P,> from level i to the continuum k is computed 
as 



ajiv) Jy 
hv 



dv [s-i] 



(A.3) 



where v, is the threshold photoionization frequency for level i, 
and the photoionization cross-sections q;,(v) for all considered 
atomic levels are taken from the TOPBASE Opacity Project on- 
line atomic database at the 'Centre de donnees astronomiques 
de Strasbourg' The photo-recombination rates P^, from the 
continuum k to level / are then given by 



. r a,(v 
^"1 — 



)By{T) 



hv 



1 - exp 



By{T) denoting the KirchhofF-Planck function at local tempera- 
ture T and frequency v. 

Cross-sections for the coUisional ionization and excitation 
by electrons are computed according to the prescriptions of 
Seaton (1962) and Van Regemorter (1962), respectively, as given 
by Allen (1976) and Cox (2000). The oscillator strengths ftj 
from Table A.2 are needed to calculate the colUsional excitation 
cross-sections. CoUisional ionization by neutral hydrogen via the 
charge transfer reaction H(li) + Li(nZ) — > Li^(li^) + W , and the 
reverse process H~ + Li^(li^) — > H(li) + \A{nl), are treated ac- 
cording to Barklem, Belyaev & Asplund (2003) for the first 7 
levels. CoUisional excitation by neutral hydrogen is ignored, as 
it was found to be unimportant for thermalizing Li i by Barklem, 
Belyaev & Asplund (2003). 

The departure coefficients bi{x,y,z) are finally obtained by 
solving the statistical equiUbrium equations locally at each grid 
point (.<:,>', z). The profile-averaged radiation field at the line 
transitions, J{vij), determines the radiative excitation rates 



Rij=BijJ{vij) 



«^line i j 



and the radiative de-excitation rates 



Rji - — Bji 

gi 



2hvh 



+ Jiyij)\ [s-'] , 



(A.5) 



(A.6) 



http : // cdsweb . u- strasbg . £r/topbase/topbase . html 



where the statistical weights of lower and upper level, g, and 
gj, and the Einstein coefficients Bij are taken from Table A.2. 
For the Une profile <py we assume a purely Gaussian distribution 
(identical for absorption and emission) with a line width that 
corresponds to the local thermal Doppler velocity plus a micro- 
turbulence of 1.5 km/s. The Doppler shift caused by the Une-of- 
sight component of the hydrodynamical velocity field is ignored 
in the present computation. For all transitions except for the res- 
onance Une, Jivij) was replaced by the mean continuum inten- 
sity Jy. For the resonance Une (transition # 1), the line opacity 
was taken into account for the computation of /(v,y), assuming 
an /-value of 0.48975 (about 2/3 of ffie value given in Table A.2 
to represent only the main component of the doublet), together 
with a typical lithium abundance of A(Li)=2.2. 

Since /(v, j) is a non-local quantity that depends in turn on 
the bi, a A-iteration is employed to obtain a consistent solution. 
Fortunately, the A iteration converges very rapidly, because even 
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the Li I resonance line affects the radiation field only marginally, 
while all the other lines are very weak. Typically, three iterations 
are sufficient to achieve convergence. 

Finally, the departure coefficients bi are used in the line for- 
mation code LinforBD to compute the non-LTE line opacity 

^low - ^up exp{-^} 

^ii„e(non - LTE) = , ' ^iir,e(LTE) , (A.7) 

l-exp{-^} 

and the Une source function 



5y,une(non-LTE) = b, 



up 



l-exp{-#} 
b\ov/-bnp expj-^} 



(A.8) 



where foiow and b^p are the departure coefficients of lower and 
upper level, respectively, and By is again the Kirchhoff-Planck 
function. 



Appendix B: Analytical fit to 3D NLTE Li abundance 
as a function of stellar parameters and EW 

B.1. Introduction 

It is often convenient to have tabular data on an irregular grid 
y = yix) condensed into a single function (even with many fit 
parameters), instead of having to interpolate between the origi- 
nal data points. Choosing a fitting function that is linear in all 
parameters (e.g., a sum of polynomials as A + Bx + Cx^ or 
sine/cosine functions) results in a comparatively simple hnear 
problem for finding these parameters, usually with one unique 
solution. However, the result might not have the desired extrap- 
olation behavior, might be "too wiggly" between the data points, 
or produce no good fit at aU. 

A more effective approach might then be to start with a prop- 
erly crafted function, possibly non-linear, and having fitting pa- 
rameters "within" the terms as in A exp(B + Cx). In this case, 
an iterative process is necessary to determine these parameters. 
This iteration requires starting values, which will affect the con- 
vergence of the process and also the solution if one is found: in 
this case, the scheme might diverge or converge to different local 
minima, depending on the starting point. However, this last issue 
is hardly a problem for our application - we are content with a 
good fit, and do not necessarily need to achieve the closest one. 
We continue, however, to have the problem of constructing the 
fitting function and choosing the initial values for the parame- 
ters. 

B.2. Method 

We used a scheme for recursive term substitution, written in IDL, 
that attempted not only to find optimum fitting parameters for a 
given function, but also to automatically determine the optimum 
functional form itself. The latter was chosen from a set of can- 
didates given by a list of terms and some construction rules. The 
usual starting point was a constant "function" with one parame- 
ter, say A, the best "fit" for which is the mean of y. In each recur- 
sion step, all parameters were successively substituted with all 
of the terms from a list (with possible restrictions), for instance 

A -> Ao +Aix, 



A -> Ao -H Al exp(A2x). 



The initial values of (Aq, Ai, A2) were derived from A by e.g., 
assuming Aq = A, Ai = 0.1 x A, A2 = 0. Applying a modified 
version of the IDL function curvefit (requiring also partial 
derivatives) provided the optimum parameter set for the given 
function and start values. 

This process was repeated with a number of random varia- 
tions to approach an optimum solution for the given function. 
The final solution was the globally best fit among all candidate 
functions and parameters. 

The main control parameters are the list of candidate terms 
and the recursion depth. For the current application, we search a 
4D function 



y = yixo, xu X2, X3) 

where ixo,xi,X2,X3) = ( log(EW), logrefr, log g, [Fe/H]), or, 
for the inverse problem, ( A(Li), logTeff , log g, [Fe/H]). We start 
with 



Ao + logio[l-exp(-10(^''-^'>)] 



to determine the equivalent width as a function of abundance, 
based on a simple absorption model for a box-shaped Une pro- 
file. In the following recursion steps, we replace the parameters 
with linear terms 



A —> Ao + Al * xi + A2X2 + A3X3, 

which finally results in higher-order polynomials instead of the 
initial coefficients. In this step, the variation in the COG with the 
stellar parameters is taken into account. In the process of apply- 
ing the term substitution, we scanned through several hundred 
possible functional forms. We decided the most suitable one on 
the basis of simplicity, and its ability to represent the numerical 
data with sufficient accuracy. 

As described in Sect. 6.1, synthetic grids were computed for 
three cases, 3D NLTE, ID LHD NLTE, and ID LHD LTE. This in 
principle leads to 6 fitting functions being computed, when both 
the EW— »A(Li) and A(Li)— >EW forms are required. We decided 
to force all three functions in either sense to have the same form, 
differing only by their parameters. We thus first searched an op- 
timum function for one case (A(Li)— >1D LHD NLTE turned out 
to be the optimum choice) and then varied only the parameters 
(and not the functional form) for the other two cases. The final 
inversion of the function was performed by hand. 



B.3. Best-fitting functions 

Best-fitting functional forms were produced both to derive A(Li) 
from EW and to derive EW from A(Li), the latter being use- 
ful, e.g., for the preparation of observations. Analytical fit for- 
mulae for both cases are presented in Eq. B.l; fitting coeffi- 
cients are listed in Table B.l. Three sets of coefficients are listed. 
The 3D NLTE set is used to fit computations where CO^BOLD 
3D hydrodynamical models are used in association with a 3D 
NLTE time-dependent spectrosynthesis. The ID NLTE and LTE 
cases refer instead to computations performed by using LHD one- 
dimensional models, with spectrosynthesis performed with or 
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Table B.l. Coefficients for the analytical fit of EW ^ A(Li) and A(Li) 





3D NLTE 


ID NLTE 


ID LTE 


Ao 


2.1744416E-I-00 


2.201 1840E+00 


2.1694977E-I-00 


Al 


3.9685178E-I-02 


-2.3708574E-I-02 


-3.1329205E-I-02 


A2 


-2.1920459E+02 


9.9754471E+01 


1.5692142E+02 


A3 


-9.8448749E+02 


1.5211893E+02 


1.9027467E+03 


A4 


2.4222436E+00 


-7.0910416E+00 


8.0236683E+00 


A 

A5 


J.U4/UoiUb+UI 


-9.jyj2UJ4b+UU 


- 1 .yjo4 j4orl+UI 


A6 


7.4822784E+02 


-1.0479124E+02 


-1.4933164E+03 


A7 


-6.0743892E-01 


1.8554645E+00 


-2.1075335E+00 


As 


-1.8920995E+02 


2.5357561E+01 


3.908771 lE+02 


A9 


-7.9977289E-02 


-7.5381881E-01 


-2.7466467E-01 


Aio 


-3.2394665E-01 


3.8146901E-01 


-4.9764225E-01 


A„ 


1.5911137E-I-01 


-2.1918788E+00 


-3.4126404E+01 


A12 


2.3078753E-02 


1.9988796E-01 


7.3497586E-02 


Al3 


8.3344564E-02 


-9.9645615E-02 


1.3085939E-01 



without the inclusion of NLTE effects, based on the same physics 
and model atoms as for the 3D NLTE case: 

gi = As + A 1 1 log Teif + A 1 2 log ^ + A 13 [Fe/H] 

22 = Afi + Qi log Teff + Ay log^ + A,^,[F c I H] 

23 = A2 + A5 log Jeff + Q2 log 5 + Aj[FelH] 
Qa=Ai + Qs log Teff + A3 log^ + A4[Fe/H] 

(B.l) 

A(LO = Al + a log Teff + A3 logg + AdFe/H] 
-l-log[-ln(l - WIO"-^")] 

logEW = Ao + log[ |1 + expC-lO^'^'^-fi")! ]• 

Units are expected to be K for temperature, cm s"^ for log 
g, and mA for EW. As for any fit, extrapolation reUabiUty is 
difficult to assess. Owing to the high computational cost of 3D- 
atmosphere model calculations, we limited our grid to the cur- 
rently available models, which obliged us to extrapolate towards 
high gravities by 0.4 dex and towards low metalhcities by about 
0.6 dex. Both extrapolations should be quite safe, since the 670.8 
nm Li doublet should be quite insensitive to both parameters. 
Because of the vanishing line opacity, atmosphere models are 
scarcely sensitive to metallicity variations below [Fe/H] = -3. 
Those who might wish to employ the presented formulae to de- 
rive Li abundances are nevertheless advised to use caution when 
extrapolating, especially in T^ff and EW. In particular, the sat- 
uration part of the COG is almost unsampled, and one cannot 
expect the present fit to reproduce it properly. Fitting to EW sig- 
nificantly above 100 pm is therefore not advisable. 



